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The  diatomic  oxides  of  Groups  IVb  and  Vb  have  been 
isolated  by  co-condensation  of  a molecular  beam  of  the 
monoxides  with  rare  gases  at  14 °K.  Zero  field  absorbance 
(ZFA)  and  magnetic  circular  dichroism  (MCD)  spectra  were 
obtained  for  all  six  molecules.  A detailed  analysis  of  the 
results  obtained  for  TiO,  ^ZrO  and  HfO  is  presented. 

An  apparatus  for  the  preparation  of  rare  gas  matrices 
containing  high- temperature  metal  oxides  has  been  construc- 
ted. ZFA  and  MCD  spectra  were  obtained  using  a computer 
controlled  system,  and  data  were  stored  on  magnetic  tape 
cassettes.  The  ZFA  spectra  were  obtained  using  a double 
beam  technique  employing  a single  detector,  two  light 
beams  modulated  at  different  frequencies,  and  two  lock-in 
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amplifiers.  The  MCD  spectra  were  obtained  using  a low  field 
electromagnet  (5500  gauss) , a photoelastic  modulator,  and 
a feedback  circuit  which  controlled  the  detector  gain  to 
correct  for  variations  in  the  background  light  intensity. 

The  ground  state  of  TiO  is  At  the  low  temperatures 

employed  in  this  study,  only  transitions  from  the  multi- 

plet  are  observed.  Since  the  multiplet  is  magnetically 

non-degenerate,  no  Cq  terms  are  seen  in  the  MCD  spectra.  The 
results  of  this  study  confirm  the  assignment  of  the  y,  y', 
a and  e systems  of  TiO  to  B^n,  C^A  and  E^n,  respec- 

tively. Several  ultraviolet  transitions  of  perturbed  high 
energy  bands  are  attributed  to  and  states.  Transitions 
attributable  to  Ti  atoms  are  observed,  and  several  assign- 
ments are  revised  in  concordance  with  the  MCD  data. 

The  ground  state  of  ZrO  is  confirmed  as  by  the  lack 

of  MCD  Cq  terms  observed  for  ZrO  bands.  The  assignment  of 
the  A system  as  and  the  C system  as  is  also 

confirmed.  The  D system  assignment  is  revised  from  ^11  X^E^ 

14-1  + 

to  E X E , and  weak  transitions  observed  at  4855  A and 
4632  A are  assigned  to  a ^ X^E"^  system.  A forbidden 
intermultiplicity  transition  is  confirmed  as  ^ X^E'*'. 

Seventeen  Zr  atomic  transitions  are  observed  in  MCD  and  ZFA 
spectra  along  with  new  ultraviolet  bands  of  ZrO  assigned  to 
high  energy  -f-  X^E^  transitions. 

The  ground  state  of  HfO  is  confirmed  as  ^E  as  well  as 

the  F and  G systems  which  are  attributed  to  ^E-^^E  or 
3 1 

E transitions  in  the  past.  The  MCD  spectra  favor  the 
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former  designation.  Similarly,  the  B system  is  either 
11  3 1 

n-^  Z or  however,  the  MCD  data  are  inconclusive  for 

this  state.  The  D and  E systems  of  HfO  are  confirmed  as 
^n-f-^Z  transitions,  and  an  ultraviolet  transition  near  2600  A 
is  proposed  as  ^n-*-^Z.  This  band  has  been  previously 
reported  in  gas  emission  spectra  but  has  not  been  seen 
previously  in  matrix  studies. 
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CHAPTER  1 


INTRODUCTION 

On  September  13,  1845,  Michael  Faraday  wrote  in  his 
journal,  ' Today  worked  with  lines  of  magnetic  force,  passing 
them  across  different  bodies  (transparent  in  different 
directions)  and  at  the  same  time  passing  a polarized  ray  of 
light  through  them."^  The  phenomenon  he  observed  was  to 
later  bear  his  name.  When  plane-polarized  light  is  passed 
through  any  transparent  substance  collinearly  with  an  exter- 
nal magnetic  field,  the  plane  of  polarization  is  rotated — 
the  Faraday  effect. 

Linearly  polarized  light  consists  of  equal  amounts  of 
right  (RCP)  and  left  (LCP)  circularly  polarized  components 
(see  Chapter  2 for  a further  discussion  of  circularly 
polarized  light) . If  the  refractive  indices  for  RCP  and  LCP 
light  are  unequal,  optical  rotation  occurs.  This  inequality 
can  be  generated  in  two  ways.  Natural  optical  activity  is  a 
result  of  low  symmetry  in  a molecule  or  in  the  unit  cell  in 
the  case  of  crystals,  and  its  wavelength  dependence  is  known 
as  optical  rotatory  dispersion  (ORD) . If  the  absorption 
coefficients  for  an  electronic  transition  are  different  for 
RCP  and  LCP,  this  is  called  circular  dichroism  (CD).  In 
natural  optical  activity  the  molecular  dissymmetry  causes  the 
electrons  to  move  in  a helical  path  on  excitation,  and 
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the  difference  in  absorption  coefficients  depends  upon  the 
handedness  of  the  helix.  The  second  way  to  bring  about  the 
inequality  is  by  an  external  magnetic  field,  H,  in  the 
propagation  direction  of  the  light.  In  magnetic  optical 
rotation  (MOR)  the  right-  and  left-handed  circular  motions 
about  H no  longer  interact  equivalently  with  the  medium,  and 
in  regions  of  absorption  the  absorption  coefficients  of  RCP 
and  LCP  light  differ  in  the  presence  of  a magnetic  field, 
giving  rise  to  magnetic  circular  dichroism  (MCD) . 

Natural  optical  activity  had  early  triumphs  in  inorganic 

and  organic  chemistry,  but  this  was  not  the  case  with  MCD 

and  MOR.  The  early  work  was  restricted  almost  entirely  to 

MOR  measurements  and  predominantly  to  transparent  spectral 

regions,  since  polarizabilities  were  most  easily  calculated 

outside  regions  of  absorption.  The  quantum  mechanical  bases 

3 

for  these  studies  were  developed  around  1930  by  Rosenfeld, 
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Kramers,  and  Serber.  The  familiar  A-B-C  terminology 
follows  the  work  of  Serber  (see  Chapter  2) . But,  it  was  not 
until  nearly  120  years  after  Faraday's  work,  that  the  modern 
methods  of  analyzing  MCD  were  developed  in  the  thesis  work 
of  Phillip  J.  Stephens  in  1964.^  The  basic  effects  causing 
magnetic  optical  activity  were  long  understood,  but  the 
important  features  of  his  work  included  improvement  in  the 
methods  of  analysis  in  regions  of  absorption  and  an  increased 
understanding  of  the  useful  information  contained  therein. 

MCD  has  now  completely  replaced  MOR  since  MCD  eliminates  many 
problems  due  to  background  rotations  from  other  electronic 
transitions  of  the  system  from  windows,  solvents,  etc.. 
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which  are  often  much  larger  than  the  effect  of  interest. 

MOR  is  contributed  to  by  all  transitions  of  the  systems  and 
therefore  is  generally  less  informative. 

MCD  can  provide  data  on  both  the  ground  and  excited 
states  of  a molecule.  In  the  case  of  the  ground  state 
properties,  there  is  some  overlap  with  magnetic  susceptibil- 
ity and  magnetic  resonance  experiments,  although  in  some 
cases  MCD  possesses  advantages  over  these  techniques,  vide 
infra.  The  prime  utilization  of  MCD  has  been  in  obtaining 
excited  state  information  which  is  unavailable  otherwise. 
Excited  state  magnetic  moments  and  g values  have  been 
obtained  in  addition  to  information  on  spin-orbit  coupling 

7 

and  Jahn-Teller  interactions. 

The  field  of  MCD  expanded  very  rapidly  during  the 

0 

1960 's,  and  a review  by  Stephens  in  1974  lists  nearly  400 

references  for  the  period  1960-1973,  whereas  a similar 
9 

review  of  the  period  from  Faraday's  time  until  1965  had 
only  133  references. 

The  Faraday  effect  is  closely  related  to  the  Zeeman 
effect;  in  fact,  MCD  is  the  difference  between  LCP  and  RCP 
Zeeman  spectra  and  provides  no  new  information  when  the 
Zeeman  components  of  a transition  are  resolved.  Hov;ever, 
the  resolution  of  Zeeman  components  is  often  difficult, 
requiring  very  high  magnetic  fields  and,  in  the  case  of 
broad  bands,  is  not  feasible  for  the  current  state  of  magnet 
technology.  The  advantage  of  MCD  is  that  it  is  a signed 
technique  and  does  not  require  high  fields;  the  present 
study  used  a very  small  electromagnet  with  a comparatively 
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weak  field  of  ~5500  gauss.  Sometimes  MCD  can  be  observed 
when  the  absorption  band  of  a sample  is  too  weak  to  be 
observed.  Typical  values  for  the  minimum  detectable  absor- 
bance is  10  ^ whereas  the  minimum  detectable  MCD  is  10“^. 

This  situation  has  been  observed  for  sharp  lines  on  para- 
magnetic species  at  low  temperatures.  An  example  occurs 
in  spin-forbidden  transitions  of  octahedral  and  tetrahedral 
Co(II)  species. 

Analytical  applications  of  MCD  have  also  been  made  for 
the  detection  and  quantification  of  various  species, 
especially  biological  molecules.  In  fact,  since  porphyrins 
exhibit  large  visible  MCD  effects,  the  technique  was  used  in 
an  attempt  to  detect  their  presence  in  lunar  rock  samples. 

Transitions  can  sometimes  be  resolved  in  MCD  and  not  in 
absorption.  A transition  overlapping  a stronger  transition 
is  resolvable  if  the  MCD  is  much  larger  for  one  of  the  tran- 
sitions due  to  the  selections  rules  for  MCD  or  if  the  transi- 
tions have  different  MCD  signs.  For  example,  tryptophan  and 
tyrosine  can  be  distinguished  from  each  other  even  though 
their  absorption  bands  overlap.  The  requirement  of  making 

the  Kramers-Kronig  transformation  of  MOR  data  to  MCD  data  in 
order  to  analyze  the  results  is  eliminated  in  MCD.^^~^^ 

The  technique  of  MCD  lends  itself  well  to  use  in 
conjunction  with  the  method  of  matrix  isolation  spectroscopy. 
Matrix  isolation  was  first  employed  by  George  Pimentel^^  in 
1954  for  trapping  unstable  reactive  intermediates  and 
involves  the  codeposition  of  a large  amount  of  solvent  gas 
with  a small  amount  of  solute  gas  species  at  low  temperature . 
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The  solvent  gas  freezes  and  traps  the  solute,  thus  isolating 
the  solute  molecules  from  each  other.  Solvent  to  solute 
ratios  are  typically  > 1000:1.  This  method  then  allows  one 
to  study  the  spectra  of  the  solute  at  leisure  and  has  a 
number  of  primary  advantages: 


1.  The  matrices  are  usually  deposited 
and  observed  at  temperatures  <25°K. 
This  simplifies  the  electronic  spectra 
in  that  "hot  bands"  are  usually  elimi- 
nated . 

2.  Except  for  a few  molecules  {H2O,  HCl , 
NH3  and  NH2)  the  rotational  structure 
of  the  spectra  is  eliminated  since  the 
species  are  no  longer  free  to  rotate 
in  the  crystal  lattice. 

3.  The  electronic  absorption  bands  are 
often  narrow  enough  to  enable  use  of 
isotopic  substitution  to  look  for  0-0 
transitions . 

4 . Reactive  species  can  be  prevented  from 
reacting  by  the  matrix,  and  higher 
concentrations  of  them  are  then 
available  for  study. 


The  basic  assumption  of  this  technique  is  that  the 
matrix  does  not  interact  with  the  solute,  and  essentially 
gas-like  information  is  obtained.  Although  this  is  a good 
first  approximation,  matrix  interactions  do  occur  which 
shift  the  energies  of  electronic  transitions,  and  the  solute 
can  be  trapped  in  multiple  sites  (different  micro-environ- 
ments) in  the  crystalline  lattice  which  leads  to  changes  in 
the  spectra.  Often  these  site  effects  can  be  detected  by 
warming  the  matrix  slightly  to  allow  controlled  diffusion 
to  occur  and  allowing  the  solute  molecules  to  settle  in  the 
more  stable  sites  thus  simplifying  the  spectra.  Lattice 
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modes  and  Jahn-Teller  distortions  due  to  interactions  with 

the  matrix  have  also  been  observed.  Recently,  Mowery  et 
, 20 

al . have  observed  Jahn-Teller  distortions  in  the  MCD  of 

matrix  of  Mg,  An  excellent,  although  somewhat  dated, 

discussion  of  matrix  isolation  is  contained  in  Beat  Meyer's 

2 1 

book  on  low  temperature  spectroscopy. 

Although  there  has  been  considerable  work  done  on 

matrices  using  optical  and  magnetic  spectroscopy,  few 

experiments  have  been  done  using  polarized  light  or  Zeeman 
22 

studies.  The  first  report  of  MCD  being  applied  to  a matrix 

. 2 ^ 

isolated  molecule  was  in  1973  by  A.  J.  Thomson,  and  only 

a handful  of  molecules  have  been  studied  since  then:  Hg,^^ 

24,25  26  20  ^ , . . 27  ^ 28  29 

C>2'  ^^2'  Xenon  halides,  benzene,  OsO^ , 

. • 30 

and  acrolein.  MCD  possesses  a special  advantage  for  study 
of  the  ground  state  magnetic  properties  of  matrix  isolated 
species.  At  present  if  the  electron  spin  resonance  spectrum 
is  undetectable,  the  magnetic  properties  are  unknown  since 
measurement  of  magnetic  susceptibilities  of  matrices  is  not 
possible  due  to  the  small  qunatities  of  material  involved. 

The  study  reported  here  involves  the  matrix  isolation 
and  MCD  of  a group  of  molecules  known  as  the  high-tempera- 
ture diatomic  oxides.  The  designation  "high-temperature" 
stems  from  the  fact  that  these  molecules  are  generated  via 
heating  the  triatomic  and  higher  oxides  to  temperatures  above 
1500°K.  The  molecules  studied  are  the  diatomic  oxides  of 
Groups  IVb  and  Vb  of  the  Periodic  table  and  are  among  the  most 
stable  molecules  known.  Their  open  shell  structure  provides 
some  interesting  visible  spectroscopy,  and  they  often  have 
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magnetic  moments  in  their  ground  and  excited  states  which 
are  amenable  to  MCD  studies. 

The  presence  of  these  oxides  has  been  observed  in  the 
stellar  spectra  of  class  M and  S stars  whose  surface  tempera- 
tures are  less  than  3000°K.^^  There  is  also  evidence  for  low 
concentrations  of  TiO,  ZrO  and  VO  in  the  comparatively  cooler 
regions  of  sunspots  and  in  X“Cygni.^^  HfO  and  TaO  have  been 
tentatively  identified  in  the  spectrum  of  R-Cygni.^^  NbO  has 
also  been  tentatively  identified  in  an  R-type  star  in  the 
constellation  Gemini.  The  spectrum  of  TiO  is  also  of 
technological  interest  since  TiO  is  observed  in  jet  engine 

exhaust  and  during  re-entry  of  space  capsules.  ZrO  is  also 
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seen  in  red  giants  and  S-type  stars. 

An  excellent  review  of  the  matrix  isolation  of  high- 
temperature  molecules  was  published  in  1978  by  Weltner,^^ 
and  his  extensive  work  in  this  field  provides  the  basis  upon 
which  this  work  rests. 

The  matrix  materials  were  Ar  and  Kr.  These  noble 
elements  are  ideally  suited  as  isolants  and  have  often  been 
the  substances  of  choice.  They  are  particularly  suited  for 
MCD  studies  since  their  crystal  structures  are  cubic  close 
packed  and  therefore  isotropic — a requirement  for  MCD. 

Neon  matrices  were  not  available  since  our  apparatus  could 
not  achieve  the  low  temperatures  required. 

The  remainder  of  this  thesis  shall  discuss  the  theory 
of  MCD,  the  details  of  the  apparatus  constructed  for  this 
experiment,  and  the  procedure  employed  in  obtaining  the 
results.  A detailed  description  of  the  results  for 
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Group  IVb  (TiO,  ZrO,  and  HfO)  is  presented.  Group  Vb 

results  are  treated  in  a more  detailed  discussion  in  the 

thesis  of  the  author's  compatriot  and  friend, 

3 8 

Robert  D.  Brittain.  Finally,  suggestions  for  further 
extensions  of  this  work  are  presented  in  Chapter  VII. 


CHAPTER  II 


THEORY  OF  MAGNETIC  CIRCULAR  DICHROISM 
Introduction 

In  this  chapter  the  pertinent  equations  which  describe 
magnetic  circular  dichroism  (MCD)  will  be  developed.  The 
approach  will  follow  the  series  of  articles  by 
Stephens , ^ emphasizing  the  approximations  and  limi- 
tations of  the  method,  A convenient  starting  point  is  a 
basic  description  of  circularly  polarized  light. 

Light  consists  of  electric  and  magnetic  fields  varying 
sinusoidally  with  the  direction  of  propagation.  In 

unpolarized  light  the  directions  of  the  electric  field 

->■ 

vector,  E,  and  the  magnetic  field  vector,  H,  are  random. 

In  plane  polarized  light  the  E and  H vectors  are  constant  in 
direction  but  vary  in  magnitude.  If  plane  polarized  light 
is  passed  through  a transparent  birefringent  medium  in  which 
the  speed  of  light  depends  upon  the  direction  of  polariza- 
tion of  the  light,  then  part  of  the  beam  will  be  delayed 
with  respect  to  the  other,  depending  upon  the  difference  in 
refractive  indices  and  the  thickness  of  the  medium.  If  the 
plane  polarized  beam  enters  the  birefringent  medium  at  45° 
to  two  perpendicular  axes  with  differing  indices  of  refrac- 
tion, equal  magnitude  components  of  the  projection  of  the 
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electromagnetic  vectors  will  be  phase  shifted  relative  to 
each  other.  When  this  phase  shift  is  equal  to  the  wavelength, 
A,  the  emergent  beam  will  be  plane  polarized  again.  If  the 
phase  shift  is  not  equal  to  A,  the  resultant  E and  H will 
trace  helical  paths  on  emergence.  Whenever  the  tracing  is 
clockwise  while  looking  toward  the  source,  the  beam  is 
considered  right  polarized  and  vice  versa  for  left  polarized. 
If  the  phase  shift  equals  the  cross  section  of  the  helix 
will  be  a circle,  and  the  light  is  referred  to  as  circularly 
polarized.  A phase  shift  different  from  an  integral 
multiple  of  causes  an  elliptical  tracing  and  the  light  is 
elliptically  polarized.  In  circularly  polarized  light  the  E 
and  H vectors  are  constant  in  magnitude  but  vary  in  direc- 
tion, but  in  elliptically  polarized  light  they  change  both 
in  magnitude  and  direction. 

A medium  exhibits  circular  dichroism  when  it  propagates 
right  and  left  circularly  polarized  light  with  unequal 
absorption.  Natural  circular  dichroism  is  a relatively  rare 
phenomenon,  whereas  magnetic  circular  dichroism  (MCD)  is  a 
property  of  all  matter.  MCD  occurs  when  a magnetic  field  is 
applied  along  the  propagation  direction  of  the  circularly 
polarized  light. 


Basic  Equations 

The  solutions  of  Maxwell's  equations  for  circularly 
polarized  light  of  angular  frequency  m traveling  in  the 
positive  z direction  in  a Cartesian  coordinate  system  yield 
for  the  electric  field  vector 
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E^(z,t)  = E°exp[io)  (t-n^z/c)  ] (i±i  j)  (1) 

and  for  the  magnetic  field  vector 

H^(z,t)  = (fi^/0)E°exp[iw  (t-fj^z/c)  ] ( + ii+])  (2) 

where  + refers  to  right  circularly  polarized  light  (RCP)  and 
- to  left  circularly  polarized  light  (LCP) , E°  is  the 
electric  field  at  z=0;  t is  the  time;  i=/^;  and  c is  the 
speed  of  light.  The  vectors  i and  j are  the  unit  vectors  in 
the  X and  y directions,  respectively.  The  complex  magnetic 
permeability  in  the  xy  plane  which  must  be  isotropic  is 
represented  by  p.  The  complex  refractive  index  is  fj^  and 
equals  p^-ik^  where  n+  is  the  real  part  of  the  refractive 
index  and  is  the  absorption  coefficient.  n is  a function 
of  the  external  magnetic  field  via  the  effect  of  the  field 
upon  k^. 

The  intensity  of  the  light,  I(z),  at  a point  z is  the 
energy  passing  through  a unit  area  normal  to  the  z axis  in 
unit  time  and  is  given  by  the  time  average  of  the  Poynting 
vector  (c/4tt)  (ExH) 

I+(z)  = (c/4tt)  Re{fi^/0}  (E°)  ^exp  (-2o)k^z/c)  (3) 

Making  the  approximation  that  the  magnetic  permeability  is 
unity  and  substituting  m = 2itv  = E/lii  where  E is  the  energy  of 
the  light  and  <h  is  Planck's  constant  divided  by  2tt,  yields 

I^(z)  = (c/4tt)  (E°)  ^exp  (-2Ek^z/^c)  (4) 

Setting  the  magnetic  permeability  equal  to  unity  excludes 
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magnetic  materials  from  consideration.  Equation  4 can  be 
further  simplified  to 

I^(z)  = I ^ (0)  exp  (-2Ek^z/l^c)  (5) 


whence 


= [Vn  + E|E^(z)  I ^][--~^(z)  ] (6) 

where  |e^(z) | is  the  magnitude  of  the  electric  field  at  z. 
The  derivative  portion  of  this  equation  is  the  energy- 
absorbed  per  unit  time  per  unit  volume  at  z and  is  related 
to  the  probability  of  absorption  per  unit  time  at  z 

at  the  energy  of  the  transition  from  state  a to  state  j , 


||(Z)  = 

0 z a,  J 


N,  P 


(z)  E 


(7) 


where  is  the  number  of  absorption  centers  of  state  a per 
unit  volume. 

By  limiting  our  discussion  to  electric  dipole  transi- 
tions and  ignoring  electric  quadrupole  and  magnetic  dipole 
interactions,  we  can  relate  the  above  probabilities  of 
absorption  to  the  electric  dipole  transition  matrix  elements 
via  time-dependent  perturbation  theory,  thereby  assuming 
that  the  effective  Hamiltonian,  Hq , is  a sum  of  independent 
components , i . e . , 


H 


o 


H"  + H, 


(8) 


where 


H°  = the  Hamiltonian  for  the  system  in  the  absence 
of  light.  The  superscript  "o"  indicates  the 
absence  of  the  radiation  field. 
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and 


Hq  = the  Hamiltonian  for  the  perturbation  due  to  the 
radiation  field. 


The  subscript  zero  indicates  the  absence  of  an  external 


magnetic  field.  This  treatment  leads  to 


P 

a->j 


( z) 


^|E^(z)  |^|<a|m^|  j>|^6  (E.^-E) 


where 


a> 


j> 


eigenstate  of 
eigenstate  of 


H°  of  energy  E . 

3. 

Hq  of  energy  E^ . 


(9) 


E . = E . - E . 

: a : a 

a = Lorentz  effective  field  correction  factor  and 
relates  the  electric  field  due  to  the  light  at 
the  absorption  center  to  the  macroscopic  E 
field . 

6 = Kronecker  delta  operator.  This  term  arises 

from  the  neglect  of  transition  lifetimes  and 
will  be  eliminated  later. 

m^  = m^jlimy  and  is  the  transition  dipole  moment 
operator  for  RCP  and  LCP  transitions. 


The  equation  for  k then  becomes 


k = ^ ^.N  |<a|m|j>|^6(E,  -E) 

± a, 3 a'  i+iJi 


(10) 


This  expression  can  be  transformed  into  more  familiar  spec- 
troscopic notation  in  terms  of  absorbance 

exp{-2Ekz/Ac)  = 10~^  - (11) 

where 


A = absorbance. 


I(z)  = the  intensity  at  point  z. 
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I (o)  = the  intensity  at  z = 0 

e = the  molar  extinction  coefficient, 
c = concentration  of  absorbers. 


This  leads  to  an  equation  for  the  absorbance 


A. 


c 

CZ  = Y4.  { ^ • 

[Na] 

E 

'±  a,j 

N 

<a|m^|  j>p6(E,  -E)  }cz  (12) 
3 a 


where 


Y+  = 


Nir^a^  logj^Qe 

250Kcti^ 


and  M is  Avogadro ' s number  and  N /N  is  the  fractional 

d. 

population  of  absorbers  in  state  a.  This  equation  assumes 
that  Beer's  law  is  obeyed  and  does  not  allow  for  non-linear 
effects  that  would  occur  at  high  light  intensities.  The 
electronic  excitations  producing  a band  are  localized,  and 
the  concentration  of  the  species  is  sufficiently  low  that 
interaction  is  negligible. 

The  dichroism  of  a system  is  the  difference  in  absor- 
bance between  LCP  and  RCP  light 


AA  = A_-A_|_ 


(13) 


Substituting  for  A_  and  A from  Equation  12  yields 


AA  r T.  . 

r “ 


N. 


[ I <a|m_| j> 


IN 

- !<a|m^|j>|2]6(E.^-E)}t: 


(14) 


where  Y=Y^=Y_  and  n^=n_.  This  equation  is  general  and 
applies  to  natural  circular  dichroism  as  observed  in 
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optically  active  molecules.  Magnetic  optical  activity  and 
MCD  are  not  perturbations  of  natural  optical  activity  but 
arise  from  different  polarizabilities  which  are  zero  when  an 
external  magnetic  field  is  absent.  Natural  and  magnetic 
optical  rotations  and  dichroism  are  additive. 

If  a system  has  no  natural  circular  dichroism  and  no 
magnetic  field  is  present,  then  AA  = 0 and  |<a|m_^|j|  = 

l<alm_|j>|  . In  this  case  we  can  define  the  zero  field 
absorbance  (ZFA)  as 

A°  = A°  = A°  = + A°)  (15) 

At  this  point  it  is  convenient  to  apply  the 
Born-Oppenheimer  approximation  (BO)  and  the  Franck-Condon 
approximation  (FC)  to  a description  of  the  ZFA.  The 
Born-Oppenheimer  approximation  assumes  that  the  zero-field 
eigenstates  can  be  written  as  a product  of  wavef unctions 

Ia  a>  = "I  (r,R)x  (R)  a = 1 to  d^ 

‘a  A A 

a 

|j,j>  = '1'^  (r,R)x.(R)  A = 1 to  d (16) 

A J ^ J u 

where  1'  indicates  an  electronic  wave  function  with  a para- 
metric dependence  upon  the  instantaneous  nuclear  configura- 
tion; X (R)  X-  (R)  sre  vibrational  wavef  unctions  and 

s J 

dependent  upon  the  nuclear  coordinates  R and  the  particular 

electornic  configuration;  r is  the  electron  coordinate; 

and  the  degeneracies  of  states  a and  j are  given  by  d^  and 

d , respectively.  The  eigenfunction  equations  are 
J 
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H A a>  - E Aa> 
o ' a a ' a 


= E.|jp>  . 


(17) 


The  Franck-Condon  approximation  is  now  employed  to  simplify 
the  transition  matrix  elements 


<A^a|m,|j,j>  = <a|j><A  |m,^|j,° 
a'i'A-'  a'±'A 


(18) 


The  lower  case  letters  refer  to  the  vibrational  wavefunc- 
tions,  and  the  superscript  zero  indicates  that  the  electronic 
matrix  element  is  evaluated  at  R = , the  equilibrium 
configuration  (the  ground  state  vibrational  functions  x 

3. 

have  significant  amplitude  only  for  R = R^  if  nuclear  motion 
is  small) . The  Zeeman  splittings  of  the  ground  and  excited 
states  will  be  independent  of  the  vibrational  level  and 
equal  to  the  pure  electronic  splittings  at  R^ . If  the  tran- 
sition is  only  weakly  allowed,  the  FC  approximation  is 
inadequate,  and  more  account  must  be  taken  of  the  R depen- 
dence of  the  electronic  wavef unction . Equation  14  now 
becomes : 


A" 

r 


= y{—  ^ 

a, A 

r\ 

Z N. 


<A  m ^ 
a ' ± 


(E  . -E)  }cz  . 


(19) 


I 

Y 


ntegrating  over 
is  a constant 


all  the  vibrational  bands  and  assuming  that 
for  all  vibronic  transitions  yields 


<A 

a 


o 

Ja" 


2 

cz. 


(20) 
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If  we  then  let 


(21) 


and 


(22) 


where 


J(E)dE  = 1 


then  the  zero  field  absorbance  is  given  by 


A° 


(23) 


Note  that  from  this  equation  we  conclude  that  the  zero  field 
absorbance  dispersion  depends  only  on  the  ground  and  excited 
state  vibrational  functions,  and  the  total  integrated 
intensity  is  dependent  on  the  electronic  wavef unctions  at 
the  equilibrium  position,  . The  integrated  intensity  is 
independent  of  the  temperature,  but  the  shape  is  dependent 
upon  the  temperature  from  the  Boltzmann  term  (N^/N)  in 
Equation  22. 

In  the  presence  of  a magnetic  field,  the  absorption 
coefficients,  k,,  exhibit  field  dependence,  and  the  eigen- 
function and  eigenvalues  of  the  system's  Hamiltonian  must  be 
obtained  as  explicit  functions  of  the  field.  This  is 
feasible  in  analytic  form  only  when  the  field  dependent 
terms  may  be  treated  by  perturbation  theory  which  assumes 
that  the  magnetic  field  energy  is  small  compared  with  the 
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separation  of  the  zero-field  states.  To  the  Hamiltonian, 

//q,  of  Equation  8 is  now  added  the  perturbation  Hamiltonian, 
H'  due  to  the  external  magnetic  field 

H ^ + H'  . (24) 

Since  nuclear  magnetons  are  three  orders  of  magnitude 
lower  than  electronic  magnetons,  we  can  ignore  the  effect  of 
nuclear  moments  and  consider  electronic  contributions  to  the 
first  order  magnetic  field  perturbation. 


H'  = + 2S  )H 

i zmc  z . z 

IX 


(25) 


where  i sums  over  all  electrons  with  mass  m and  charge  e. 

H is  the  magnitude  of  the  magnetic  field  directed  along  the 

z axis;  is  the  projection  on  the  z axis  of  the  angular 

i 

momentum  of  the  ith  electron  and  similarly  for  the  projection 

of  the  spin  S . Summing  over  all  electrons  leads  to 
i 


H'  = -M  HEB (L  + 2S  )H 
z z z 


(26) 


where  \x ^ is  the  electronic  magnetic  moment  along  z,  and  3 is 
the  electronic  Bohr  magneton  = 4.6681  x 10~^  cm~^/gauss.  L 
and  S^  are  the  corresponding  total  orbital  and  spin  angular 
momenta  along  the  z axis.  This  equation  is  recognized  as 
describing  the  Zeeman  splittings  induced  by  the  magnetic 
field  and  leads  to  population  changes  in  the  ground  state. 
The  Boltzmann  term  in  Equation  22  can  now  be  described  in 
terms  of  the  magnetic  moments  and  the  external  magnetic 
field : 
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N d^exp{-<A  |l  +2S  |a  >6H/kT} 
g _ g g ' z z ' g 

N “ Zd  exp{-<A  |L  +2S  |A  >3H/kT} 
g g g I z z i g 


where  a sums  over  the  Zeeman  states  of  A.  For  high  tempera- 
tures and/or  low  field  strengths  the  exponentials  can  be 
expanded  via 


- X 

e 


1 - X + 


2 

X 

yr 


(28) 


whereby 


^ = Sr‘i-"\l^^  + 2SjA_^>6H/kT)  . (29) 

If  one  assumes  that  the  Zeeman  shift  merely  changes  the 
position  of  the  absorption  band  along  the  energy  axis  and 
does  not  change  the  shape  of  the  band,  the  shifted  lineshape 
function  can  be  expressed  in  terms  of  the  unshifted  function 
by  a Taylor's  series  expansion.  This  approach  is  called  the 
rigid  shift  approximation  and  is  especially  valid  for  broad 
bands  where  the  Zeeman  shift  is  small  compared  with  the  band 
width.  The  assumption  that  the  band  moves  rigidly  will 
usually  be  a good  zero-order  approximation,  but  may  require 
modification  since  little  is  known  about  the  variation  of 
line  shapes  in  magnetic  fields  and  we  are  implicitly  ignoring 
any  magnetic  field  perturbations  to  the  transition  moment 
operators . 

Finally,  the  manner  in  which  the  circularly  polarized 
absorption  is  modified  by  the  magnetic  field  using 
Equations  13,  19,  27  and  the  above  considerations  leads  to 
an  expression  for  the  MCD: 
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A A = +y{A. 


6fl 

( C 1 

651 

+ 

f } (3H)  cz 


where 


(30) 


The  three  terms  in  this  equation  are  referred  to  as  the 
8q  and  terms  and  the  shape  of  the  MCD  can  be 
deconvolved  into  these  component  parts.  Conversely,  calcu- 
lation of  these  terms  can  be  used  to  predict  the  sign  and 
shape  of  the  MCD  dispersion.  The  maximum  contributions  to  the 
MCD  of  the  three  terms  are  related  as 


A 


1 


:6, 


_Z  Z 
AW‘  kT 


(31) 


where  T is  the  bandwidth  of  an  electronic  transition; 


AW  is 
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the  order  of  magnitude  of  an  electronic  energy  gap;  k is  the 
Boltzmann  constant;  and  T is  the  absolute  temperature.  Z is 
the  Zeeman  energy. 

In  examining  Equation  30  one  notes  that  the  existence  of 
a Cq  term  indicates  that  the  ground  state  is  degenerate,  and 
the  Cq  term  depends  upon  the  magnetic  moment  of  the  ground 
state  and  the  circularly  polarized  transition  moments. 

Since  the  Cq  term  is  divided  by  kT,  for  Cq  0 the  temperature 
dependence  of  the  MCD  depends  upon  the  magnitude  of  the 
term.  Note  also  that  the  term  is  multiplied  by  the 
lineshape  function,  f. 

The  presence  of  an  term  indicates  that  the  ground 
state  and/or  the  excited  state  is  degenerate  and,  since  it 
is  multiplied  by  the  derivative  of  the  lineshape  function, 
depends  strongly  on  the  bandwidth  of  the  transition.  Thus 
the  MCD  experiment  provides  information  on  the  magnetic 
moments,  symmetries,  and  electronic  wavef unctions  of  the 
ground  and  excited  states. 

The  6^  term  contribution  to  the  MCD  depends  upon  mixing 
of  the  electronic  wavef unctions  via  the  magnetic  field 
perturbation  and  is  the  chief  cause  of  the  Faraday  effect 
being  present  in  all  substances,  since  there  will  always  be 
some  mixing,  however  slight,  with  other  states  k. 

The  subscripts  on  the  three  MCD  terms  in  Equation  30 
and  in  Equation  21  refer  to  the  moments  of  the  MCD  and 
ZFA  curves  that  yield  these  values.  In  1965  Henry, 
Schnatterly  and  Slichter  introduced  the  moments  method  of 
analyzing  MCD  of  color  centers. 


The  moments  of  a spectrum 
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are  averaged  properties  where  the  nth  moment  of  the  ZFA  and 
MCD  are  defined  as 

rO  /'  A° 

<A°/  = J -f(E-E°)'^dE  (32) 

and 

<Aa/°=  j ~(E-E°)'^dE  (33) 

where  E°  is  the  energy  of  the  center  of  the  band  about  which 
the  moment  is  taken.  The  terms  are  related  to  the  moments  by 

<A°>o  = <e°>cz  = 326.6P^cz  (34) 

<AA>o  = <Ac°>cz  = 326.6[8q  +^J  (3H)cz  (35) 

<AA>^  = <Ae>^cz  = 326 .6A^  (3H)  cz  (36) 

Thus  one  can  obtain  the  values  for  the  terms  from  the 
experimental  curves  and  compare  these  with  the  calculated 
terms  from  Equation  30.  In  order  to  define  the  moments 
properly,  the  band  must  be  completely  resolved  from  neigh- 
boring bands  at  all  fields.  The  moments  must  be  taken  over 
the  entire  set  of  vibronic  transitions  allowed  in  conjunction 
with  the  electronic  excitation. 

For  overlapping  systems  the  MCD  remains  an  additive 
function  of  the  individual  contributions  and  allows  one  to 
separate  bands  that  overlap  in  absorption  spectra  but  have 
different  signs  or  magnitude  of  the  component  MCD  terms. 

The  effect  is  observed  for  TaO  in  these  studies  where 
oppositely  signed  terms  allow  separation  of  previously 
unresolved  band  systems.  If  the  two  overlapping  bands  are 
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very  close  in  energy  and  the  magnitude  of  the  Cq  terms  are 
comparable,  then  the  MCD  will  have  a "pseudo-A"  term  with 
the  appearance  of  a derivative  when  in  actuality  it  is 
merely  an  overlap  phenomenon  (see  Figure  1) . A similar 
"pseudo-A"  term  can  sometimes  be  observed  when  strongly 
interacting  states  lead  to  oppositely  signed  Bq  terms. 

Selection  Rules  for  MCD  and  Angular 
Momentum  of  Diatomic  Molecules 

Photons  possess  well-defined  values  of  angular  momentum 
and  are  emitted  and  absorbed  by  systems  which  are  character- 
ized by  well-defined  values  of  their  angular  momenta. 

Since  angular  momentum  is  conserved,  the  selection  rules 
characterizing  transitions  involving  photons  with  angular 
momenta  must  take  this  into  account. 

Let  us  begin  by  considering  the  transitions  occurring 
in  an  isolated  atom  with  total  angular  momentum  J.  When  no 
magnetic  field  is  present  the  state  is  degenerate  to  the 
extent  of  2J+1  corresponding  to  the  number  of  distinct 
eigenvalues  m^ . The  field  about  the  atom  is  spherically 
symmetric,  and  there  is  no  preferential  direction  for  the 
angular  momentum.  When  an  external  magnetic  field  is  intro- 
duced, this  degeneracy  is  lifted  via  the  space  quantization 
imposed  by  the  field  direction.  Figure  2 shows  the  Zeeman 
splitting  of  the  eigenstates  of  atoms  with  total  angular 
momentum  J=0  and  J=1  in  a magnetic  field  and  the  three 
possible  transitions  allowed. 

The  probability  for  the  absorption  of  a photon  of 
electric  dipole  radiation  for  a transition  from  states  a to  j 
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Pseudo-A^  terms  generated  by  oppositely  signed 
Cq  or  Bq  terms  for  overlapping  transitions. 


Figure  1. 


Figure  2.  Magnetic  field  splittings  and  polarized 
transitions  to  an  atomic  P-state. 


Transition  a is  RCP . 

Transition  b is  polarized  parallel  to  the 
magnetic  field. 


Transition  c is  LCP. 
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is  proportional  to  the  square  of  the  transition  dipole 
moment  matrix  element  along  the  direction  of  the  photon 
polarization. 

(u*m.  ) E (u*<a|m|j>)  (37) 

J a 

->-  , 

where  u xs  a unit  vector  determining  the  polarization  of  the 

light,  and  m^^  is  the  transition  electric  dipole  moment.  We 

can  now  derive  the  selection  rules  for  the  three  transitions 

of  Figure  2 by  considering  the  three  possible  polarizations 

of  light  (x,y,z)  on  Equation  37  and  the  properties  of  the 

4 2 

wavef unctions  as  described  by  the  spherical  harmonics. 

First  considering  the  polarization  parallel  to  the 
field,  the  z direction.  Equation  37  vanishes  for  all  cases 
except  where  AJ  = ±1  and  Am  = 0.  This  situation  is  of  no 

u 

interest  in  MCD  since  the  z polarization  is  the  light 
propagation  direction  in  MCD  and  is  unavailable. 

When  considering  the  two  possible  polarizations  of  the 
light,  u^  and  u^ , when  the  light  is  propagating  along  z,  it 
is  more  convenient  to  examine  the  two  linear  combinations 
u^±iu^  corresponding  to  the  two  possible  circular  polariza- 
tions. The  + combination  is  right  circular  polarized,  and 
the  - is  left  circular  polarized.  For  RCP , Equation  37 
vanishes  except  for  Am  = m -m  = -1  and  conversely,  the 
selection  rule  for  LCP  is  Am  = +1,  and  AJ  = ±1  for  both 

U 

cases.  From  these  considerations  we  can  conclude  that  tran- 
sition b in  Figure  2 will  not  be  observed,  whereas  transition 
a will  correspond  to  RCP  and  transition  c to  LCP.  Thus  in 
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this  case  circularly  polarized  light  corresponds  to  electric 
dipole  photons  having  J=1  and  ra  = +1  for  LCP  and  m = -1 
for  RCP.  Although  these  selection  rules  were  derived  for 
the  case  of  atoms,  the  law  of  conservation  of  angular  momen- 
tum makes  the  absorption  by  any  system  subject  to  similar 
selection  rules. 

V7e  shall  now  turn  to  the  special  considerations  that 
apply  to  the  angular  momenta  of  diatomic  molecules.  In 

—y  , 

atoms  the  electronic  orbital  angular  momentum,  L,  is  a 
constant  of  the  motion  of  the  electrons  due  to  the  spherical 
symmetry  of  the  force  field.  In  a quantum  system  which  does 
not  have  spherical  symmetry,  the  orbital  angular  momentum  is 
not  an  integral  of  the  motion.  In  diatomic  molecules  the 
force  field  is  axially  symmetric  about  the  internuclear 
axis,  and  the  component  of  the  orbital  angular  along  this 
axis  is  a constant  of  the  motion.  L processes  in  a classical 
sense  about  the  axial  electric  field  of  the  nuclei  with 
constant  component  m - L,  L-l,...-L.  Unlike  motion  in  a 
magnetic  field,  reversing  the  motion  of  the  electrons  about 
the  electric  field  does  not  change  the  energy,  and  m and 

1j 

-m  are  equivalent.  Therefore,  it  is  more  appropriate  to 

L 

classify  the  electronic  states  of  diatomic  molecules 

according  to  the  absolute  value  of  m , and  the  quantum 

L 

number  of  interest  in  diatomics  is  designated 

A = |m^  I . (38) 

The  value  of  A is  represented  in  the  term  symbol  by  E,  II,  A, 
and  (J)  according  to  A = 0,1, 2, 3,  respectively. 
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Except  for  Z states  all  diatomic  states  are  doubly 
degenerate.  External  space  quantization  can  occur  via  the 
introduction  of  an  external  magnetic  field  to  lift  this 
degeneracy,  whereby  a state  would  split  into  two  differ- 
ent m^  states  rather  than  the  three  states  corresponding  to 

the  three  values  of  m as  observed  in  atomic  states.  This 

L 

is  illustrated  in  Figure  3. 

Just  as  in  atoms,  the  spin  of  the  electrons  can  add 
vectorially  to  yield  a total  resultant  spin,  S.  In  the  case 
where  A ^ 0,  there  is  an  internal  magnetic  field  due  to  the 
orbital  motion  of  the  electrons  which  causes  space  quantiza- 
tion along  the  internuclear  axis.  The  classical  precession 
of  S about  the  internuclear  axis  has  a constant  projection 
designated  Z which  can  have  the  values 

Z = S,  S-1, . . . ,-S  . (39) 

In  contrast  to  A,  Z can  have  positive  or  negative  values. 

To  carry  the  comparison  between  atoms  and  diatomic 
molecules  further,  the  total  angular  momentum  in  atoms  is 
designated  J and  is  the  vector  sum  of  L and  S.  The  corre- 
sponding quantum  number  in  diatomic  molecules  is  simply 
the  algebraic  sum  of  A and  Z,  since  they  both  lie  on  the 
internuclear  axis,  and  is 

n = I A + Z I . (40) 

The  symbol  for  a diatomic  term  is 


Figure  3.  Magnetic  field  splittings  in  atoms  and 
diatomics  with  L or  A = 1. 
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where  2S+1  is  the  multiplicity,  and  A = E , II , A , $ for  values 
of  A = 0 , 1 , 2 , 3 . 

An  electronic  terra  with  a given  A ^ 0,  splits  into  a 
2S+1  multiplet  of  components  of  different  energies  due  to 
the  interaction  of  Z and  A,  spin-orbit  coupling.  The  elec- 
tronic energy  of  a multiplet  term  is  given  to  a first 
approximation  by 

T = T + A(AE)  (41) 

e o 

where  T^  is  the  term  value  when  spin  is  neglected.  A is  a 
constant  for  a given  multiplet,  increases  rapidly  with  the 
number  of  electrons  as  in  atoms,  and  may  be  either  positive 
or  negative.  Figure  4 shows  the  multiplet  splittings  and 
vectors  for  the  ground  state  of  TiO.  As  long  as  A ^ 0, 
each  of  the  multiplet  components  is  doubly  degenerate,  even 
those  with  n = 0 and  will  split  in  a magnetic  field  as  in 
Figure  3 except  in  those  situations  where  the  values  of  the 
orbital  angular  momentum  and  the  spin  angular  momentum 
cancel  each  other  to  give  no  net  magnetic  moment.  An  example 

3 

of  this  cancellation  occurs  in  the  A^  state  of  TiO,  vide 
infra . 

When  the  interaction  of  the  spin  and  orbital  angular 
momenta  are  so  strong  that  A and  Y.  are  no  longer  valid 
quantum  numbers,  the  selection  rules  presented  above  must 
be  modified.  This  strong  coupling  often  occurs  in  species 
containing  nuclei  with  high  atomic  numbers,  and  Q.  is  the 
only  valid  quantum  number  in  these  species.  The  MOD 
selection  rule  then  becomes  Am^^  = 


±1  and  is  known  as 
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Figure  4 
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Spin-orbit  coupling  in  diatomic  molecules. 
Relative  orientations  of  A,  Z and  are  shown. 
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Hund's  case  c.  When  the  coupling  is  weak,  as  in  species 
containing  light  nuclei,  the  MCD  selection  rules  are  = ±1 

and  Am^  = 0.  This  later  situation  is  called  Hund's  case  a. 
These  two  cases  are  the  extreme  situations,  and  often  mole- 
cules will  display  properties  intermediate  between  case  a 
and  case  c. 

To  illustrate  how  the  various  magnetic  field  splittings 
give  rise  to  the  MCD  and  Cq  terms.  Figure  5 shows  the 
MCD  transitions  for  the  cases  where  either  the  ground  state 
or  excited  state  is  split.  The  term  arises  when  the 
excited  state  or  ground  state  is  split  by  the  magnetic 
field.  A positive  A^  term  is  shown  where  the  positive-going 
lobe  is  to  the  higher  energy  corresponding  to  the  absorption 
of  LCP  light  being  more  energetic  than  RCP  light.  Two  cases 
are  illustrated  for  the  case  where  the  ground  state  is 
split.  When  the  temperature  is  high  the  populations  of  all 
the  ground  substates  will  be  equal  and  the  MCD  dispersion 
will  be  as  illustrated  in  the  middle  spectrum.  As  the 
temperature  is  lowered,  the  population  in  the  lowest  mag- 
netic sublevels  begins  to  increase  according  to  Boltzmann 
distribution  statistics.  The  lobe  for  the  higher  sublevels 
decreases  while  the  lobe  for  the  lower  levels  increases, 
the  resulting  dispersion  is  illustrated  by  the  dashed  curve 
in  the  spectrum  on  the  right. 
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CHAPTER  III 


MATRIX  ISOLATION  APPARATUS 

This  section  contains  a detailed  description  of  the 
apparatus  constructed  for  matrix  isolation  including  the 
electronics  and  optical  systems  used  in  the  absorption  and 
MCD  experiments.  Figure  6 shows  a schematic  representation 
of  the  furnace  and  matrix  isolation  device.  A rather  compact 
furnace  was  required,  because  it  had  to  fit  within  the  con- 
fines of  a small  Alpha  Model  4600  4-inch  adjustable  gap 
electromagnet.  The  distance  between  the  Knudsen  cell  and 
the  cold  window  was  3.5  inches,  which  caused  some  problems 
due  to  the  radiation  flux  from  the  Knudsen  cell,  vide  infra. 

Originally,  a large  furnace  design  was  used  which  had 
the  bulk  of  the  furnace  outside  of  the  magnet,  and  the  dis- 
tance between  the  Knudsen  cell  and  the  cold  window  was 
approximately  10  inches.  This  design  proved  unsatisfactory, 
because  this  longer  distance  allowed  a greater  dispersion  of 
the  molecular  beam  effusing  from  the  Knudsen  cell  with 
concommitantly  weaker  absorbing  matrices. 

The  molecular  beam  was  generated  in  a resistively 
heated  Knudsen  cell  made  of  either  tungsten  or  tantalum. 

These  materials  were  suitable  for  high  temperature  studies 
because  of  their  refractory  nature;  however,  some  difficul- 
ties did  exist.  The  tantalum  cells  were  made  from  lengths 
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of  0.25  inch  outer  diameter  tubing  with  a wall  thickness  of 
0.035  inch.  At  the  high  temperatures  required  for  formation 
of  the  diatomic  oxides  studied,  Ta  rapidly  reacted  with  the 
oxides,  and  often  the  cell  failed  via  cracking  before  the 
required  temperatures  could  be  maintained.  Tantalum  cells 
were  used  only  in  the  TaO  experiments . 

The  more  expensive  and  fragile  tungsten  cells  were  used 
with  the  other  experiments  because  of  the  inertness  of 
tungsten.  These  cells  were  expensive--$30  per  cell  from 
Chemetal.  Tungsten  was  also  very  brittle  and  care  had  to  be 
taken  when  putting  the  effusion  hole  in  the  cell.  This  hole 
(-0.05  inch)  had  to  be  ground,  not  drilled;  and  even  with 
care  to  avoid  thermally  shocking  the  cell,  they  seldom 
lasted  more  than  three  runs.  The  only  end  cap  material 
available  for  the  Knudsen  cells  was  tantalum,  which  proved 
to  be  convenient  since  a small  amount  of  TaO  observed  in  the 
matrix  was  often  used  as  an  internal  standard  for  the  sign 
of  the  MOD. 

The  Knudsen  cell  was  mounted  by  tantalum  straps  to  two 
water  cooled  copper  electrodes.  One  of  these  electrodes  was 
electrically  isolated  from  the  furnace  flange,  and  the  other 
was  electrically  connected  to  the  entire  outside  of  the 
furnace;  therefore,  one  needed  to  be  cautious  when  using 
the  magnetic  shutter  that  the  magnet  did  not  cause  arcing 
across  the  electrodes.  Up  to  300  amp,  60  Hz,  alternating 
current  was  passed  through  the  small  area  of  the  cell  wall 
enabling  operation  at  2500°  C.  This  high  current  was 
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R 

obtained  by  means  of  a Variac  variable  transformer  and  a 
step-down  transformer. 

Cooling  of  the  furnace  was  achieved  by  water  flowing 
through  0.25  inch  copper  tubing  soldered  to  the  outer  shell 
of  the  furnace.  During  depositions  at  temperatures  above 
2000°  C,  it  was  necessary  to  wrap  the  furnace  with  wet 
cotton  strips  and  pour  liquid  nitrogen  on  the  furnace  in 
order  to  maintain  the  integrity  of  the  furnace  "0"  ring 
seals.  The  temperature  of  the  cell  surface  was  measured 
using  a Leeds-and-Northrup  optical  pyrometer  with  subsequent 
surface  emissivity  correction  for  the  type  of  cell  material. 
The  magnetic  shutter  was  closed  except  during  temperature 
measurements  to  prevent  the  formation  of  mirror-like  films 
on  the  viewport  which  affected  the  temperature  readings. 

A large  proportion  of  the  heating  of  the  outside  of  the 
furnace  was  due  to  radiation  from  the  Knudsen  cell.  Very 
little  heat  transfer  was  from  conduction  through  the  furnace 
since  the  pressure  was  maintained  below  4 x 10  torr.  In 
order  to  prevent  this  radiant  flux  from  overloading  the 
cooling  capacity  of  the  cold  window  (4  watts  at  14°K) , it 
was  necessary  to  place  a water  cooled  heat  shield  between 
the  Knudsen  cell  and  the  cold  window.  This  heat  shield  was 
made  by  sandwiching  a spiral  of  0.25  inch  copper  tubing 
between  two  0.125  inch  copper  plates.  A 0.125  inch  hole 
was  drilled  through  both  plates  which  also  served  to  improve 
the  collimation  of  the  molecular  beam,  thus  reducing  the 
thermal  load  of  molecules  hitting  parts  of  the  cold  window 
assembly  other  than  the  actual  sample  window.  Careful 


42 


alignment  of  the  heat  shield  hole  and  the  Knudsen  cell  hole 
was  necessary  to  ensure  that  the  molecular  beam  impinged 
correctly  on  the  cold  window. 

The  molecular  beam  was  cocondensed  with  an  inert  gas  on 

the  surface  of  the  cold  window.  The  temperature  of  the 

R 

window  was  kept  at  -14°  K using  an  Air  Product  Displex 
Model  CS  202  closed  cycle  helium  refrigerator.  The  temper- 
ature of  the  cold  window  was  monitored  using  a Kp/Au-0.07  at 
% Fe  thermocouple  mounted  near  the  bottom  of  the  cell  window 
assembly.  This  assembly  was  made  from  a high  purity, 
oxygen-free  copper  block  with  a h"-28  threaded  stud  for 
mounting  to  the  second  stage  of  the  refrigerator.  Indium 
gaskets  were  used  to  give  good  thermal  conduction  between 
the  window  holder  and  the  refrigerator.  These  gaskets  were 
also  used  between  the  cold  window  and  its  holder  and  between 
the  heat  shield  and  the  first  stage  of  the  refrigerator. 

The  cold  window  was  a disk  of  CaF2  1 inch  in  diameter 
and  0.25  inch  thick.  A Suprasil  II  quartz  window  was 
originally  used;  however,  this  material  was  unsatisfactory 
for  MCD  since  it  exhibited  a strongly  temperature-dependent 
birefringence,  probably  induced  during  cooling  by  strain  due 
to  differences  in  the  coefficients  of  thermal  expansion  of 
copper  and  quartz.  Although  CaF2  could  exhibit  strain  bire- 
f ringence--it  was  used  as  a material  for  photoelastic  modula- 
tors--perhaps  the  strain  of  cooling  was  more  easily  relieved 
since  CaF2  was  somewhat  softer  than  quartz.  The  windows  on 
the  outside  of  the  vacuum  shroud  were  2-inch  diameter 
Suprasil  I quartz,  and  no  birefringence  could  be  detected 
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due  to  their  presence  in  the  light  beam.  The  temperature  of 
the  cold  window  was  controlled  by  a 10  watt  variable  duty 
cycle  heater  in  a feedback  loop  with  the  thermocouple 
signal.  There  was  a 2°  K temperature  gradient  between  the 
top  and  the  bottom  of  the  window  holder.  The  minimum  temper- 
ature obtainable  was  10°  K. 

The  inlet  nozzle  for  the  isolant  gas  was  made  from  a 
16  gauge  Yale  stainless  steel  needle  which  was  silver- 
soldered  to  a 0.25  inch  diameter  stainless  steel  tube.  This 
tubing  was  silver-soldered  through  a hole  in  the  wall  of  a 
short  length  of  brass  tubing  which  joined  the  furnace  and 
the  vacuum  shroud  of  the  cryostat.  A 0.25  inch  vacuum  quick 
connect  joined  the  inlet  nozzle  assembly  to  a 0.25  inch 
stainless  steel  tube  leading  to  a glass  inert  gas  manifold. 
The  inert  gas  was  fed  into  the  vacuum  system  via  a Nupro 
Extra  Fine  Metering  Valve.  The  gas  was  passed  through  a 
liquid  nitrogen  cooled,  stainless  steel  trap  as  an  extra 
precaution  to  remove  any  contaminants.  The  entire  gas 
manifold  was  diffusion  pumped  to  1 x 10~^  torr  and  baked  at 
60°  C under  vacuum  for  several  days  before  the  high  purity 
gas  was  introduced  into  the  three  liter  reservoirs  via  a 
Hoke  0.25  inch  stainless  steel  needle  valve. 

A diagram  of  the  rare  gas  manifold  is  shown  in 
Figure  7.  Note  that  there  are  two  reservoir  bulbs  enabling 
rapid  change  between  different  isolant  gases.  The  diffusion 
pump  was  a 3 inch  Varian  M2  loaded  with  100  ml  of  Dow  Corning 
704  silicone  oil.  Ace  Glass  Company  Teflon^  valves  were 
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used  for  all  other  connections.  The  flow  rate  of  the  inert 
gas  was  measured  using  a standard  mercury  side-arm  manometer. 

During  deposition  of  a matrix,  the  pressure  in  the 
system  had  to  be  low  enough  that  the  mean  free  path  of  the 
effusing  molecules  was  longer  than  the  distance  from  the 
furnace  to  the  cold  window.  The  pressure  was  maintained 

_ 5 

below  4 X 10  torr  by  a diffusion  pumping  system  shown  in 
Figure  8.  The  lowest  pressure  obtained  with  this  system  was 

- 7 

4 X 10  torr  when  the  furnace  was  not  in  operation.  The 
pressure  was  monitored  by  an  ionization  gauge  and  thermo- 
couple gauge  connected  to  a Granville-Phillips  Series  270 
Gauge  Controller.  The  ionization  gauge  was  mounted  as  close 
as  possible  to  the  Knudsen  cell.  The  thermocouple  gauge  was 
mounted  on  the  0.25  inch  line  used  for  rough  pumping  of  the 
system  before  the  diffusion  pump  was  opened  to  the  system 
via  a 2-inch  diameter  Consolidated  Vacuum  Company  gate 
valve.  A stainless  steel  liquid  nitrogen  trap  was  between 
the  gate  valve  and  the  2 inch  Consolidated  Vacuum  Company 
diffusion  pump,  charged  with  80  ml  of  Dow  Corning  704  sili- 
cone oil.  The  system  was  rough  pumped  to  a pressure  below 
200  microns  before  the  gate  valve  was  opened.  Both  of  the 
roughing  pumps  were  Welsh  Scientific  Company  Duo  Seal 
Model  1400. 


Spectroscopic  Apparatus 


The  cryostat  was  rotated  so  that  the  cold  window  faced 
the  furnace  and  nozzle  during  deposition.  When  spectroscopic 
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experiments  were  performed,  the  window  was  rotated  90°  to 
intercept  the  light  beam.  Both  MCD  and  absorption  experi- 
ments were  performed  with  approximately  15  minutes  required 
for  interchange  between  the  two  modes  of  operation. 

This  section  deals  with  the  optics  and  electronics  for 
the  absorption  and  MCD  experiments.  Finally,  a description 
of  the  computer  system  and  interface  is  presented. 

A diagram  of  the  absorption  optics  is  given  in 
Figure  9.  The  apparatus  was  designed  to  operate  as  a 
"pseudo-double  beam"  spectrometer.  The  reference  beam  was 
not  passed  through  a matched  cell,  but  since  the  various 
optical  components  and  the  matrix  gases  were  transparent 
over  the  spectral  regions  studied,  this  presented  little 
problem.  Scattering  from  the  matrices  did  cause  some 
background  drift  at  the  shorter  wavelengths,  but  the  lamp 
emission  spectrum  was  effectively  nulled  over  most  of  the 
spectral  range. 

Spectra  were  obtainable  over  the  range  of  2500  to 

O 

8500  A.  There  was  a strong  emission  peak  in  the  spectrum  of 

O 

the  high  pressure  xenon  lamp  around  8400  A which  exceeded 
the  dynamic  range  of  the  lock-in  amplifiers  used.  This 
region  gave  poor  nulling  of  the  background  during  absorption 
experiments.  The  lamp  was  a 300  watt  Eimac^  made  by  Varian 
Associates.  This  lamp  operated  at  a pressure  of  115 
atmospheres  of  xenon  and  had  a built-in  parabolic  reflector 
and  a sapphire  window  yielding  a well-collimated  intense 
beam.  The  lamp  was  housed  in  a small  hood  for  removal  of 
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Figure  9.  Absorption  experiment  optics. 
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the  ozone  generated  during  operation.  This  hood  was  attached 
to  a Spex  0.75  meter  Czerny-Turner  Spectrometer  with  a 
1200  line/mm  grating  blazed  at  3000  A with  an  f-number  of 
6.8.  The  slits  were  set  for  a spectral  band  pass  of  1.5  to 
3 A which  is  about  the  limit  of  resolvable  structure  in  Ar 
and  Kr  matrices. 

The  beam  emerging  from  the  monochromator  was  passed 
through  a 3 inch  quartz  lens  (f/3) . This  lens  spread  the 
beam  so  that  it  could  pass  through  two  sets  of  slots  on  a 
spinning  wheel.  The  wheel  was  driven  by  a hysteresis- 
synchronous  motor  at  1800  rpm.  The  wheel  contained  5 outer 
slots  and  9 inner  slots  giving  chopping  frequencies  of 
150  Hz  and  270  Hz.  These  frequencies  were  selected  to  mini- 
mize ac  powerline  pick-up  since  they  were  not  integral 
multiples  of  the  ac  line  frequency. 

The  light  beam  chopped  at  270  Hz  was  passed  through  the 
matrix  and  onto  the  photocathode  of  an  EMI  9683  QB  photomul- 
tiplier tube  in  a Fact-50^  Peltier  effect  cooler.  The  spec- 
tral response  of  the  tube  was  extended  S-20,  and  the  tempera- 
ture of  the  tube  could  be  maintained  at  -30°  C for  minimiza- 
tion of  dark  current.  The  light  beam  chopped  at  150  Hz  was 
reflected  around  the  matrix  by  two  Edmund  Scientific  alumi- 
nized front-surface  mirrors  and  onto  the  same  photomultiplier 
tube . 

The  chopping  wheel  assembly  was  equipped  with  two  small 
1.5  volt  light  bulbs  and  two  phototransistors  located  on 
opposite  sides  of  the  chopping  wheel  from  the  lights  such 
that  the  signal  from  the  phototransistors  were  used  as 
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reference  signals  for  the  respective  chopping  frequencies. 
These  signals  were  amplified  by  two  RCA  3140  operational 
amplifiers.  Both  the  light  bulbs  and  the  amplifiers  were 
powered  by  batteries.  A diagram  of  the  absorption  elec- 
tronics is  given  in  Figure  10. 

The  current  output  from  the  photomultiplier  tube 

_ 7 

(~10  amp)  was  converted  to  a voltage  by  a 15  resistor 
connected  to  the  ground.  Part  of  the  signal  was  amplified 
by  a factor  of  100  using  an  Analog  Devices  52  k low-drift 
operational  amplifier  and  sent  to  an  oscilloscope  for  moni- 
toring the  wave  shape  during  the  experiment.  The  rest  of 
the  signal  was  presented  in  parallel  to  the  inputs  of  two 
lock-in  amplifiers:  an  Ithaco  Model  353  Phase-lock  Amplifier 

and  an  Ithaco  Dynatrac  391A  Lock-in  Amplifier.  The  Model  353 
amplifier  was  locked  to  the  150  Hz  reference  frequency  from 
the  chopping  wheel,  and  its  output  was  proportional  to  the 
lamp  emission  spectrum,  monochromator  dispersion  character- 
istics, and  the  photomultiplier  spectral  response.  The 
Model  391A  was  locked  to  the  270  Hz  reference  signal,  and 
its  output  was  proportional  to  the  factors  mentioned  above 
plus  the  absorption  spectrum  of  the  matrix.  Both  outputs 
were  connected  to  a Semiconductor  Circuits,  Inc.,  Log 
Amplifier  which  took  the  logarithm  of  both  signals  and 
performed  an  analog  subtraction  to  give  an  absorbance  output. 
This  output  was  connected  in  parallel  to  the  computer 
interface,  vide  infra,  and  the  input  of  a Texas  Intruments, 
Inc.,  Servo/Riter  II  Potentiometric  Recorder. 
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The  MCD  experiments  were  performed  with  the  same  light 
source  and  monochromator.  The  mirror,  lens,  and  chopping 
wheel  for  the  absorption  experiments  were  replaced  by  a 
cylindrical  quartz  lens  (f/7)  to  collimate  the  beam  which 
was  passed  through  a Glan-Thomson  prism  oriented  at  an  angle 
of  45°  to  the  modulation  axis  of  a Morvue  Electronic  Systems 
PEM-3  photoelastic  modulator.  All  current  high  sensitivity 
MCD  instrumentation  was  based  upon  phase  modulation  using 
either  electro-optic  modulators  (EOM)  or  photoelastic  modu- 
lators (PEM) . All  EOM  employed  ADP/KDP  crystals  and  required 
high  voltages  applied  longitudinally.  The  angular  aperture 
was  small,  and  the  transparent  range  was  somewhat  limited. 

A PEM  was  chosen  because  of  the  larger  angular  aperture  and 
lower  voltages.  See  Figure  11  for  a diagram  of  the  MCD 
optics . 

An  Alpha  Model  4600  electromagnet  with  a 4 inch 
adjustable  pole  face  gap  was  used  in  the  MCD  experiments. 

This  magnet  had  a 0.75  inch  hole  collinear  with  the  magnetic 
field,  as  required  for  MCD,  and  was  mounted  on  tracks  of 
case-hardened  steel  via  open  ball-bushings.  This  mobile 
mounting  allowed  the  magnet  to  be  freely  moved  in  and  out 
of  position  as  required. 

Figure  12  shows  the  MCD  electronics.  The  signal  from 
the  photomultiplier  tube  was  connected  to  an  Ithaco  Model  193 
Lock-in  Amplifier  which  was  locked  to  the  50  kHz  reference 
signal  from  the  PEM.  The  bipolar  output  of  this  lock-in 
^plifier  was  level  shifted  to  yield  only  positive  voltages 
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Figure  11.  Optics  for  MCD  experiments. 
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Figure  12.  Electronics  for  MCD  experiments. 
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suitable  as  chart  recorder  input  and  input  to  the  computer 
interface.  The  output  from  the  Analog  Devices  52  k opera- 
tional amplifier  was  sent  to  an  error  amplifier-feedback 
circuit . 

This  feedback  circuit  (see  Figure  13)  was  similar  to 

4 3 

that  of  Krausz  and  Cohen,  but  employs  a Bertan  PMT-20 
(Option  3)  programmable  high  voltage  power  supply.  The 
reader  should  be  cautioned  that  there  are  errors  in  the 
schematic  in  the  Krausz  and  Cohen  article  and  the  circuit 
presented  there  will  not  function  as  drawn.  The  feedback 
circuit  controlled  the  gain  of  the  photomultiplier  tube  in 
an  inverse  relation  to  the  dc  output  of  the  photomultiplier 
tube.  The  signal  from  the  PM  tube  consists  of  a relatively 
small  50  kHz  ac  component  riding  on  a larger  dc  component 
which  is  proportional  to  the  lamp  emission  spectrum  and  the 
sample  absorbance.  Maintaining  a constant  dc  level  by  the 
feedback  circuit  automatically  corrects  for  these  background 
effects  in  the  MCD  spectrum. 

Both  the  absorption  and  MCD  experiments  were  computer 
controlled.  The  original  design  employed  a KIM  I system 
from  MOS  Technology,  Inc.,  but  this  system  was  later  replaced 
by  the  more  powerful  Commodore  PET  2001  microprocessor.  The 
scanning  of  the  wavelength  was  accomplished  by  pulses  sent 
from  the  computer  to  a SLO-SYN^  Model  MOGI-FD-301  stepping 
motor  attached  to  the  hand-scan  control  on  the  monochromator. 
This  motor  stepped  200  times  per  revolution  allowing 
0.25  A/step. 


Figure  13.  Feedback  Circuit  for  MCD 
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The  data  were  digitized  every  angstrom  using  a Datel 
Systems,  Inc.  ADC-Ek8B  analog-to-digital  converter  with  8 
bit  resolution.  Because  the  PET  computer  had  a rather 
modest  memory  capacity  (8000  bytes),  the  maximum  number  of 
data  points  that  could  be  taken  and  stored  by  the  computer 
was  2000.  The  microprocessor  was  based  on  6502  architecture 
which  required  data  storage  in  blocks  containing  no  more  than 
256  data  points;  therefore,  the  data  were  stored  in  250  point 
blocks  which  were  labeled  alphabetically  from  A to  H.  After 
collecting  a spectral  scan,  the  data  were  stored  on  magnetic 
cassette  tapes.  A complete  listing  of  the  program  for  data 
collection  and  storage  is  found  in  the  appendix. 

The  computer  also  controlled  the  photoelastic  modulator 
synchronized  with  the  wavelength  drive.  The  wavelength  of 
quarter-wave  retardation  by  the  modulator  depended  upon  the 
amplitude  of  the  stress  applied  to  the  quartz  crystal  which 
was  linearly  dependent  upon  an  input  voltage  to  the  modula- 
tor. The  program  calculated  the  required  voltage  after  each 

O 

0.25  A step,  and  this  output  was  sent  to  the  modulator  via 
a Datel  Systems,  Inc.,  DAC-IC8BC  digital-to-analog  converter. 
This  same  device  was  also  used  when  a second  program  read 
data  from  the  magnetic  cassette  tapes  for  output  to  the 
chart  recorder  or  to  the  Nicolet  1180  computer.  This 
program  is  also  listed  in  the  appendix. 

A program  was  written  in  1180  BASIC  for  the  input  of 

O 

data  fields  up  to  2000  A wide  into  the  Nicolet  1180  computer, 
storage  on  the  Diablo  disk  system,  and  scaled  plotting  on 
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the  digital  plotter  for  hard  copy  files.  The  program  dis- 
played the  spectrum  on  the  digital  oscilloscope  and  plotted 
it  using  scale  lengths  determined  by  the  operator.  This 
program  is  also  listed  in  the  appendix  along  with  a second 
program  in  1180  BASIC  used  to  assemble  a full  8500  to  2500  A 
scan  from  data  blocks  stored  on  the  disk  and  for  calibration 
of  the  final  spectra.  The  assembled  spectrum  was  then 
stored  on  the  disk  for  later  plotting  using  another  program 
which  is  also  listed  in  the  appendix. 


CHAPTER  IV 


EXPERIMENTAL  TECHNIQUES 

This  chapter  is  concerned  with  the  experimental 
procedure  employed  in  obtaining  a matrix  isolated  sample  of 
the  high  temperature  oxide  of  interest.  The  general  tech- 
niques were  very  similar  for  all  the  oxides  studied;  the 
only  real  difference  being  the  furnace  temperatures  employed. 
The  various  calibration  procedures  are  also  presented. 

As  mentioned  in  Chapter  III,  the  material  employed  for 
the  Knudsen  cells  was  either  tungsten  or  tantalum.  These 
materials  were  made  by  a reductive  extrusion  process  in  the 
presence  of  an  atmosphere  of  hydrogen,  which  became  trapped 
in  the  metal  interstices.  Therefore,  the  cells  were  out- 
gassed  by  heating  to  2200 °C  under  vacuum  conditions  for 
approximately  30  minutes.  The  release  of  hydrogen  was 
readily  monitored  by  observing  the  pressure  in  the  furnace. 

The  cell  was  then  loaded  through  the  effusion  hole  with 
the  required  metal  oxide.  The  cell  was  loaded  as  fully  as 
possible  since  some  of  the  oxide  invariably  blew  out  into 
the  vacuum  chamber  during  pump  down  and  warm-up  of  the 
furnace.  Any  dust  remaining  from  previous  experiments  was 
removed  from  the  vacuum  system  by  washing  with  Kimwipes 
soaked  with  acetone.  The  furnace  flange  and  Displex^ 
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shroud  0-rings  were  thoroughly  cleaned  and  greased  with 
Apiezon-N  . 

The  Knudsen  cell  was  then  mounted  via  tantalum  straps 
to  the  water-cooled  copper  electrodes  making  certain  the 
effusion  hole  was  aligned  with  the  holes  in  the  water-cooled 
heat  shield  and  the  CaF^  window.  This  alignment  was  visually 
checked  by  sighting  the  Knudsen  cell  from  the  CaF^  window 
while  heating  the  cell  to  a dull  red  glow.  Then  by  turning 
off  the  current  to  the  furnace,  the  effusion  hole  in  the 
cell  could  be  observed  as  a spot  glowing  more  brightly  than 
the  surrounding  cell  wall  for  several  seconds.  This  effect 
was  due  to  the  fact  that  the  interior  of  the  Knudsen  cell 
cooled  more  slowly  than  the  exterior,  and  this  technique 
helped  insure  proper  alignment.  This  method  was  only 
employed  after  the  system  was  under  vacuum  conditions. 

The  furnace  was  first  rough  pumped  for  5 to  10  minutes 
using  the  bypass  line  to  isolate  the  diffusion  pump  and 
liquid  nitrogen  trap.  Rough  pumping  was  continued  until  the 
system  pressure  was  below  200  microns  as  monitored  by  the 
thermocouple  gauge.  The  gate  valve  to  the  diffusion  pump 
was  then  gradually  opened;  the  bypass  line  was  closed;  and 
the  system  was  diffusion  pumped  below  2 x 10~^  torr.  Stan- 
dard acetone  spraying  techniques  were  employed  to  check  for 
leaks . 

The  furnace  cell  was  then  gradually  heated  to  the 
deposition  temperature  while  maintaining  the  pressure  below 
2 X 10  ^ torr.  Too  rapid  warm-up  caused  blow-out  of  material 
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from  the  cell  and  led  to  weakly  absorbing  matrices;  this 
preheating  usually  took  several  hours  to  accomplish.  During 
this  time  the  Displex  was  rotated  so  that  the  heat  shield 
protected  the  clean  CaF2  window  from  material  effusing  from 
the  furnace.  This  material  could  be  observed  as  a film  on 
the  heat  shield. 

The  CaF^  window  had  to  be  cleaned  with  acetone  between 
subsequent  experiments  in  order  to  insure  good  matrices. 
Merely  warming  the  window  to  remove  a matrix  without  also 
washing  with  acetone  was  not  sufficient  cleaning  since  the 
next  matrix  deposited  was  invariably  cloudy. 

After  preheating,  the  cell  was  often  cooled  to  800°C 
and  left  overnight  at  this  temperature.  The  cell  was  then 
heated  to  200°C  below  the  deposition  temperature;  the 
Displex  compressor  was  turned  on  to  cool  the  window  to 
14 °K.  The  rare  gas  cold  trap  was  then  filled  with  liquid 
nitrogen,  and  the  leak  valve  opened  so  that  1 to  3 mmol/hr 
flowed  into  the  system.  The  CaF2  window  was  rotated  so  that 
the  rare  gas  could  impinge  on  the  cold  surface  for  10  minutes 
before  the  furnace  temperature  was  increased  to  the  deposi- 
tion temperature. 

The  flow  rate  for  the  rare  gas  was  kept  low  to  insure 
clear  matrices  and  avoid  polarization  scrambling  in  the  MCD 
experiments.  The  ratio  of  oxide  to  rare  gas  was  estimated 
to  be  1000:1.  Deposition  continued  45  to  60  minutes,  and 
many  of  the  matrices  showed  pronounced  color.  The  furnace 
temperature  was  then  reduced  to  800 °C,  and  the  rare  gas  was 
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allowed  to  flow  for  5 minutes  longer.  The  CaF2  window  was 

then  rotated  90°  so  that  its  absorbance  could  be  measured. 

If  the  absorbance  of  the  strongest  bands  was  below  0.5, 

further  deposition  was  carried  out  until  the  absorbance  was 

in  the  range  of  0.5  to  0.95.  The  maximum  signal-to-noise 

ratio  in  MCD  experiments  occurred  when  the  absorbance  was 
44 

0.87.  If  the  absorbance  was  deemed  sufficient,  the  furnace 
was  allowed  to  cool  to  room  temperature. 

Both  MCD  and  absorption  experiments  were  performed 
before  and  after  annealing  at  30  to  33°  K for  30  minutes. 
Spectra  were  taken  at  14  and  24°  K.  If  the  spectra  exhibited 
any  temperature  dependence,  they  were  recorded  at  2 to  3°  K 
intervals  between  14  and  24°  K. 

O 

The  monochromator  was  calibrated  to  within  0.5  A using 
the  lines  from  a Spectra  Physics  4W  argon  ion  laser.  Wave- 
length scanning  was  always  from  low  to  high  energy  insuring 

. O 

maximum  phototube  sensitivity  and  avoiding  a 1 A backlash 
error  in  the  stepping  motor  assembly. 

The  MCD  spectra  were  recorded  at  5.5  kilogauss  as 
determined  by  an  F.  W.  Bell  Model  640  Gaussmeter.  Zero 
field  dichroism  spectra  were  recorded  to  determine  the  zero 
line  in  the  MCD.  With  the  magnet  power  supply  turned  off 
there  was  a residual  field  of  62  gauss. 

The  MCD  apparatus  was  calibrated  for  each  matrix  spectra 
obtained  using  aqueous  solution  of  the  dichroic  complex 
[Co (en) ^]Cl • (d- tartrate)  with  an  absorbance  of  1.00  at 

O 

4690  A as  measured  on  a Cary  14  spectrometer.  The  procedure 

4 5 

was  similar  to  that  of  Tacon,  and  the  circular  dichroism 
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(CD)  of  this  solution  was  measured  at  4930  A using  the  same 

conditions  under  which  each  MCD  was  obtained.  McCaffery 
46 

and  Mason  measured  the  CD  and  absorption  spectra  of  this 
compound.  The  absorption  band  at  4690  A has  a CD  maximum 
at  4930  A with  AA  /A  equal  to  0.0225.*  The  compound  was 

synthesized  using  the  method  of  Broomhead^^  and  had  the  same 
specific  rotation,  [a]^,  of  +103°  as  reported.  This  same 
solution  was  used  to  calibrate  the  absorbance  experiment. 
Equation  42  belov/  was  used  in  the  calibration  of  the 
digitized  MCD  data. 


AA  = 


0.0202 


f (sam) 

(sam) -I^ (0) • 

f Utd) 

(std)-I^ (0) 

(42  ) 


where 


f(sam)  = lock-in  amplifier  sensitivity  for  MCD  of 
the  sample 

f(std)  = lock-in  amplifier  sensitivity  for  CD  of 
the  standard 


I, (sam)  = digitized  value  for  sample  MCD  at  wave- 
length A 

(0)  = digitized  value  for  the  zero  field  value 
at  wavelength  A . 


* Value  reported  in  McCaffery  and  Mason's  article. 
used  in  Tacon's  Thesis  0.02C2. 


Value 


CHAPTER  V 


RESULTS  OF  TiO  EXPERIMENTS 

This  chapter  and  the  next  are  concerned  with  the 
experimental  results  obtained  for  the  group  IVb  diatomic 
oxides.  The  results  for  ZrO  and  HfO  will  be  discussed 
together.  The  spectra  of  ZrO  and  HfO  are  similar  since 
their  ground  state  configurations  are  the  same  but  differ 
from  that  of  TiO.  This  reversal  of  configuration  also 
occurs  in  group  Vb  in  which  VO  differs  from  NbO  and  TaO. 

Table  1 lists  the  ground  state  configurations  of  groups  IVb 
and  Vb. 

Although  titanium  is  relatively  abundant  in  the  earth's 
48 

crust  (0.6  %)  and  TiO  is  one  of  the  most  stable  compounds 

4 9 

known--its  dissociation  energy  is  159.8+3  kcal/mole  — the 
diatomic  oxide  of  titanium  is  seldom  encountered  on  earth. 

The  preferred  form  is  Ti02  which  occurs  in  the  minerals 
rutile,  anatase  and  brookite. 

Even  though  relatively  rare  terrestrially,  TiO  is  one 
of  the  most  widely  spectroscopically  studied  transition 
metal  oxides.  With  the  exception  of  ScO,  TiO  is  the  simplest 
transition  metal  oxide  with  d electrons  involved  in  chemical 
bonds.  Any  adequate  theory  of  chemical  bonding  must  be  able 
to  account  for  the  role  that  d electrons  play.  The  spectra 
of  relatively  cool  class  M and  S stars  are  dominated  by 
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TABLE  1 


Ground  State  Configurations  of  the 
Diatomic  Oxides  of  Groups  IVb  and  Vb 


Oxide 

Spectroscopic 

Term 

Configuration 

TiO 

'a 

R 

(6a) 

ZrO 

(a  ) 

HfO 

(0^) 

VO 

(6^a) 

NbO 

(6a^) 

TaO 

(6a^) 
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absorption  lines  attributable  to  TiO;  therefore  any  theories 
of  the  thermodynamics  of  stellar  atmospheres  must  include 
consideration  of  the  role  of  TiO.  It  is  possible  to  use  the 
rotational  profile  of  TiO  to  determine  the  temperature  of 
stellar  atmospheres.  The  visible  spectrum  of  TiO  has  been 
studied  extensively  since  the  early  work  of  Lowater^^  on 
the  Y system  published  in  1927  and  Christy's  work^^  on  the 
a system  in  1929.  Since  that  time,  well  over  60  articles 
have  appeared  in  the  literature  concerning  the  spectrum  of 
the  molecule.  Tracing  the  changes  in  state  assignments 
through  the  literature  from  1927  to  the  present  demonstrates 
the  complexity  involved  in  the  analysis  of  the  spectrum  of 
this  simple  oxide.  For  example,  the  a and  y systems'  assign- 
ments have  been  revised  several  times.  In  the  late  1920 's, 

the  a system  was  considered  to  be  a transition,  and 

3 3 S2 

the  Y was  1-  II.  In  the  early  1950 's,  Phillips  revised 

3 3 

the  assignment  of  the  y system  to  A-  n by  analysis  of 

spectra  taken  using  a grating  spectrometer  21  ft  long.  In 
53 

1954,  Uhler  suggested  that  the  ground  state  of  TiO  was 
on  the  basis  of  comparison  with  the  spectra  of  ZrO. 

The  assignments  had  to  be  revised  again,  and  the  presently 
accepted  assignment  for  the  a system  is  C A-X  A.  The  MCD 
results  presented  here  confirm  these  assignments. 

Carlson  and  Moser^^  and  Carlson  and  Nesbet  ^ carried  out 
the  first  ab  initio  calculations  on  a heteronuclear  molecule 
containing  a transition  metal.  The  TiO  molecule  was  chosen 
for  study  since  it  was  the  best  characterized  diatomic  of 
this  type  with  the  fewest  number  of  electrons.  At  that  time 
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the  cost  of  calculation  of  Hartree-Fock  wavef unctions  for 
even  the  simplest  compounds  containing  third  row  elements 
was  prohibitively  expensive;  therefore,  they  employed  a 
very  limited  basis  set  representing  the  normally  occupied 
atomic  orbitals  and  those  orbitals  of  higher  angular  momentum 
which  substantially  contributed  to  the  representation  of 
polarization  and  distortion  of  atomic  valence  orbitals. 
Nineteen  basis  functions  were  used  that  had  been  previously 
optimized  as  a minimal  self-consistent  field  basis  for  the 
atoms. These  calculations  confirmed  Uhler's  prediction 
that  the  ground  state  of  TiO  was  ^A^(6a).  This  assignment 
was  in  contrast  to  a predicted  (6^)  ground  state  proposed 
by  Berg  and  Sinanoglu  based  on  the  ligand  field  splitting 
of  the  d-orbital  degeneracies  for  the  metal  ion,  Ti^"^(3d^), 
in  the  field  of  the  oxide  ion.  The  term  was  actually 

predicted  by  Carlson's  calculations  to  be  very  high  in 
energy  due  to  large  electron  repulsion  effects  and  has  not 
been  observed  experimentally.  A simple  molecular  orbital 
diagram  for  the  ground  state  of  TiO  is  presented  in  Figure  14. 
Later  calculations  by  Carlson  and  Moser^^  indicated  that  in 
the  ground  state  configuration  of  TiO,  (6o) , the  6 

R 

electron  density  is  located  primarily  on  the  Ti  atom,  and 
the  o density  is  located  primarily  on  the  0 atom. 

The  two  lowest-lying  singlet  states,  ^A  and  ^Z'^,  were 
attributed  to  the  configurations  (6a)  and  (a^) , respectively. 
The  ^Z”^  state  was  observed  earlier  by  Peterson^^  and 

and  the  ^A  state  by  Phillips. 


Peterson  and  Lindgren 
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Figure  14.  TiO  MO  diagram. 


72 
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The  prediction  of  a ground  state  was  confirmed  by 

6 2 

Weltner  and  McLeod  in  the  first  matrix  isolation  study  of 

• X 8 

TiO  in  1965.  By  employing  0-substitution  and  comparison 

with  gas  phase  results,  they  were  able  to  establish  (0-0) 

band  positions  and  vibrational  progressions  for  the  a,  y and 

y'  systems  in  neon  and  argon  matrices. 

The  TiO  matrices  reported  in  this  dissertation  were 

prepared  in  argon  under  similar  conditions  using  Ti02  powder 

(Alfa,  99+%  optical  grade) . Knudsen  cell  temperatures  were 

in  the  range  of  2150  to  2250°K  with  argon  flow  rates 

2.5  mmol/hr.  The  a,  y and  y'  systems  were  observed  in 

addition  to  the  weak  c system  later  reported  by  McIntyre, 

6 3 

Thompson  and  Weltner.  The  ZFA  band  positions  in  argon 
are  listed  in  Table  2,  and  a representative  ZFA  spectrum 
for  the  a,  y and  y'  systems  is  presented  in  Figure  15. 

The  a bands  compose  the  shortest  wavelength  grouping 
and  clearly  demonstrate  the  potential  utility  of  MCD  in 
assignment  of  symmetry.  Since  the  multiplet  splitting  of 
the  ground  state  is  regular,  the  ground  multiplet  is 

K 

3 . 3 

A^.  The  next  higher  energy  state  is  the  A^,  which  is 

—1  3 3 —1 

96.4  cm  above  the  A^  state.  The  A^  state  is  101.1  cm 
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above  the  A^  state.  Assuming  a Boltzmann  distribution 
at  24°K,  only  3 molecules  in  1000  would  be  in  the  A^  state 
and  7 molecules  in  a million  would  be  in  the  ^A,^  state. 
Therefore,  the  only  absorptive  transitions  observed  under 
these  conditions  will  originate  from  the  A^^  state.  For 
case  a coupling,  ^A^  is  a rather  special  state  since  its 
magnetic  moment  will  be  zero  due  to  the  opposition  of  the 
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TABLE  2 

Absorption  Bands  of  TiO  in  an  Argon  Matrix 


Band 

v' 

A,  A* 

o -|- 

A,  a' 

-1 

V , cm 

Y 

0 

7253 

7254 

13785 

1 

6839 

6837 

14626 

2 

6470 

6469 

15458 

3 

6142.5 

6135 

16300 

Y' 

0 

6254 . 5 

6254 

15990 

1 

5941 

5940 

16835 

2 

5657 

5655 

17683 

3 

5400.5 

5400 

18519 

a 

0 

5241.5 

5241 

19080 

1 

5023.5 

5024 

19904 

2 

4825 

4826 

20721 

3 

4645 

4634 

21580 

4 

4478 

4480 

22321 

e 

0 

8394 

8390 

11919 

1 

7798 

7795 

12829 

p 

8166 

8160 

12255 

AG  ' 


v'+h' 


-1 

cm 


841 
832 

842 


845 

848 

836 


824 

817 

817 

741 


910 


* See  References  61  and  62 
t Experimental  data. 
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spin  and  orbital  angular  momenta,  i.e,,  <^A  |l  +2S  I > = 0 

1 z z ' 1 

A similar  situation  arises  in  states  which  also  have  zero 

magnetic  moment.  Since  TiO  contains  comparatively  light 
nuclei  with  an  atomic  number  of  22  for  Ti  and  8 for  0,  one 
would  expect  case  a coupling  which  would  cause  all  transi- 
tions ^A-^A  to  be  forbidden  in  MCD.  In  case  c coupling  only 
one  transition  would  be  allowed,  i.e.,  ^A^-* — ^A^  where 
Am^  “ Figure  16  shows  the  magnetic  field  splitting  for 

these  states  and  the  only  allowed  MCD  transition  in  case  a or 
case  c.  Figure  17  shows  the  ZFA  spectrum  recorded  in  this 
region.  The  lack  of  an  MCD  signal  is  good  evidence  for 
case  a coupling  in  TiO.  Phillips^^  also  notes  that  in  the 
gas  phase  spectra  all  of  the  triplet  systems  in  TiO 
correspond  to  Hund ' s case  a.  The  small  positive  MCD  peak 

O 

at  4835  A is  due  to  the  presence  of  Ti  atoms  which  are 
present  in  the  matrix,  vide  infra.  This  Ti  atomic  absorp- 
tion is  responsible  for  the  distortion  observed  in  the 
(2-0)  band  of  the  a system.  This  transition  was  observed  in 

ft  ft 

a matrix  isolation  study  of  Ti  atoms  by  Gruen  and  Carstens 

3 0 3 

and  was  assigned  by  them  to  the  z D^-f — a F^  transition.  This 

peak  exhibited  the  temperature  dependent  behavior  expected 

for  the  C-term  characteristic  of  a ^F^  ground  state. 

Laser  excitation  into  the  system  of  TiO  was  used  by 
6 7 

Brom  and  Broida  in  1975  to  measure  the  relative  energy 

spacing  between  the  singlet  and  triplet  manifolds.  The 

value  obtained  was  3500  cm  ^ between  the  isoconf igurational 
13 

a A and  X A states  which  was  in  good  agreement  with  theoret- 
ical predictions  but  was  in  conflict  with  the  1952  results  of 
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Figure  16.  Magnetic  field  splittings  and  MCD  allowed 
transitions  in  the  a bands  of  TiO. 
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6 8 

Phillips.  The  earlier  work  was  based  on  population  anal-- 
yses  of  gas  phase  spectra  and  stated  that  an  energy 
difference  ranging  from  274  to  804  cm~^  with  an  average 
difference  of  581  cm  ^ existed  between  the  low-lying  a^A 
level  and  the  ground  state  which  was  considered  to  be  n 
at  that  time. 

Several  unsuccessful  attempts  were  made  to  use  excita- 
tion with  a Spectra  Physics  4W  Argon  ion  laser  to  observe 
magnetic  circularly  polarized  emission.  The  small  bore  of 
the  magnet  and  the  long  distance  required  between  the  sample 
and  the  photomultiplier  tube  resulted  in  the  unsuccessful 
completion  of  these  experiments.  The  photomultiplier  tube 
had  to  be  at  least  20  in  from  the  matrix  because  of  the 
large  effect  of  the  magnetic  field  on  the  gain  of  the  tube. 
Other  attempts  to  observe  magnetic  circularly  polarized 

emission  using  a more  compact  superconducting  magnet  were 

69 

also  unsuccessful. 

After  the  a system  the  next  lower  energy  band  system 

observed  in  an  argon  matrix  was  the  y'  system.  This  system 

7 0 

was  fxrst  discovered  by  Coheur  in  1943  in  gas  phase 
spectra  taken  using  an  exploding  titanium  wire.  The  y' 
bands  are  often  much  stronger  on  astronomical  plates  of 
M-type  stars  than  on  laboratory  plates^^  and  are  the  most 
intense  bands  observed  in  the  matrix  spectra.  The  y' 
system  lies  in  an  extremely  crowded  region  in  the  gas  phase 
spectra,  and  the  analysis  therefore  has  been  a difficult 
task.  Many  of  the  rotational  lines  are  still  unassigned. 
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The  ZFA  spectrum  of  the  y'  system  is  illustrated  in 
Figure  15,  and  the  MCD  spectrum  is  presented  in  Figure  18. 

The  assigned  transition  for  the  y'  system  is  . 

R 

Since  spin-orbit  coupling  occurs  in  the  B^IT  state,  the 

multiplet  structure  of  the  state  must  be  considered  in  the 

calculation  of  the  MOD  terms.  Phillips^^  gives  a multiplet 

spacing  of  21.6  cm  ^ between  the  lowest  energy  ^11^  and 

multiplets  and  16.2  cm~^  between  the  and  the  highest 

3 

energy  multiplets.  The  only  transition  allowed  in 

case  a coupling  is  the  and  is  illustrated  in 

Figure  19.  The  forbidden  transition  B^II^^X^A^  has  been 
observed  in  gas  phase  spectra.  In  MCD  spectra  involving 
case  a coupling  this  transition  is  forbidden  since  Am^  ^ 0, 
and  in  case  c coupling  the  Am^  = ±1  selection  rule  is  also 
broken.  The  B^II^-^-X^A^  transition  is  also  forbidden  in  MCD 
for  case  a coupling  since  Am^  ^ 0.  The  calculation  of  the 
value  for  the  MCD  term  A^/Pq  from  equations  21  and  30  yields 
a value  of  +13  for  the  B^n^-<-X^A^ transition  which  agrees  with 
the  experimentally  observed  sign  of  the  MCD  observed  for 
this  system.  Table  3 lists  the  positions  of  the  positive 
and  negative  going  peaks  in  the  MCD  for  the  v'  = 0,1, 2, 3 of 
the  y'  system  of  TiO  in  an  argon  matrix. 

The  broad  and  sharp  structure  present  in  the  y ' system 
is  illustrated  in  an  expanded  spectrum  of  the  (0,0)  band  in 
Figure  20.  This  band  can  be  considered  as  composed  of  two 
sharp  positive  terms  and  one  broad  positive  A^^  term 
originating  from  different  sites  in  the  matrix.  The  other 
possible  interpretation  is  that  the  structures  might  be  due 
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Figure  19.  Magnetic  field  splittings  and  MCD  allowed 
transitions  for  the  y and  y'  bands  of  TiO. 
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TABLE  3 


V ' 


Peak 

Positions  of  Positive 

and  Negative 

of  the 

Y ' System  of  TiO  in 

an  Argon  Matrix 

A (A) 

1 

e 

o 

1 

Sign 

Shape 

6290 

15898 

Broad 

6264 

15964 

- 

Sharp 

6258 

15980 

+ 

Sharp 

6256 

15985 

- 

Sharp 

6253 

15992 

+ 

Sharp 

6248 

16005 

+ 

Broad 

5969 

16753 

Broad 

5948 

16812 

- 

Sharp 

5944 

16824 

+ 

Sharp 

5941 

16832 

- 

Sharp 

5937 

16844 

+ 

Sharp 

5930 

16863 

+ 

Broad 

5684 

17593 

Broad 

5663 

17658 

- 

Sharp 

5654 

17687 

- 

Sharp 

5644 

17718 

+ 

Broad 

5404 

18505 

Broad 

5442 

18376 

+ 

Broad 
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to  the  splittings  of  the  multiplets  ^II^,  and 

Although  the  11^  state  has  been  observed  in  gas  phase 
spectra  and  the  splitting  of  the  sharp  peaks  is  close  to 
the  multiplet  spacing,  the  MCD  transitions  to  the  and 

3 

11  are  forbidden  and  a calculation  of  the  MCD  A^/V  terms 
from  equations  21  and  30  yields  -16  and  -36/  respectively. 

The  interpretation  of  these  structures  as  due  to  site  effects 
on  a single  transition  is  also  supported  by  the  fact  that  the 
sharp  structures  were  not  observed  in  all  matrices,  and  the 
relative  intensities  changed  on  annealing  the  matrix. 

Different  site  effects  are  often  noted  in  matrix  isola- 
tion experiments.  The  site  splittings  noted  here  are  among 
the  sharpest  yet  observed  for  MCD  spectra.  These  unique 
energy  sites  in  the  Ar  matrix  are  probably  due  to  the  discon- 
tinuous and  highly  structured  nature  of  the  Ar  film.  When 
Ar  is  deposited  from  the  vapor  state  at  low  temperatures  the 
crystals  formed  are  of  poor  crystalline  quality  with  very 

small  grains,  which  contain  many  twins,  stacking  faults,  and 
71 

dislocations.  The  films  are  defective  because  the  number 

of  individually  nucleated  crystals  rises  dramatically  as  the 

temperature  is  lowered  at  constant  pressure.  Kovalenko  and 
72  . 

Bagrov  studied  the  diffraction  of  electrons  from  Ar  films 
deposited  at  different  temperatures.  The  electron  diffrac- 
tion may  be  interpreted  in  terms  of  small  coherently  diffrac- 
ting blocks  where  the  block  size  is  roughly  the  distance 
between  dislocations  in  the  crystal.  For  Ar  deposited  at 
16-23°K  the  block  size  was  estimated  to  be  70-80  A and 
20-30  A for  deposition  at  5-15°K. 
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Assuming  a block  size  of  30  A for  deposition  at  14  K, 
the  isolation  temperature  used  for  TiO,  it  is  possible  to 
calculate  the  relative  number  of  substitutional  sites  on  the 
surface  and  in  the  interior  of  these  blocks.  The  total 
number  of  atoms  in  a face-centered  cubic  lattice  of  Ar  30  X. 
on  an  edge  is  approximately  1000.  There  are  370  substitu- 
tional sites  on  the  outer  surface  layer  and  468  sites  within 
one  layer  of  the  surface.  Therefore,  approximately  half  of 
the  sites  are  on  the  surface  of  these  small  Ar  crystals. 
Annealing  at  30-35°K  allows  the  diffusion  of  TiO  to  different 
sites,  but  the  temperature  must  be  raised  to  48°K  for  the 
size  of  these  Ar  crystals  to  change. 

Very  broad  peaks  are  observed  to  lower  energies  of  the 
Y'  and  Y bands.  These  peaks  are  most  likely  due  to  the 
presence  of  other  species  present  in  the  matrix.  These 

73 

peaks  were  also  present  in  the  spectra  taken  by  Weltner. 

Some  of  the  pertinent  observations  are  listed  below: 

1.  The  intensity  of  the  broad  peaks  relative  of  the  y and 
Y ' bands  in  the  ZFA  and  MOD  spectra  varies  between 
matrices . The  spectrum  of  Ti  atoms  dominate  below 
4000  A,  and  the  matrices  which  have  strong  bands  for  Ti 
atoms  also  have  strong  bands  for  these  broad  features. 

2.  The  MCD  bands  of  the  broad  peaks  display  the  character- 
istic inverse  temperature  dependence  of  Cq  terms, 
implying  that  the  ground  state  is  magnetically  degener- 
ate. Figure  21  shows  the  MCD  temperature  dependency  of 
the  broad  peak  at  3574  A. 

3.  The  relative  intensity  of  the  broad  peaks  to  the  y and 
Y'  bands  varies  greatly  between  the  ZFA  and  MCD  spectra. 
In  the  MCD  spectra  the  broad  peaks  are  much  more  intense 
than  in  the  ZFA  spectra  which  would  further  support  the 

observation  since  the  magnitude  of  a term 
contribution  is  much  greater  than  an  term. 
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4.  The  sign  of  the  term  for  the  broad  peaks  associated 

with  the  y'  system  is  negative  and  positive  for  the 

Y system. 

5.  As  the  temperature  of  the  Knudsen  cell  during  deposition 
was  increased  from  2140  to  2260°K,  the  proportion  of 
these  broad  peaks  also  increased. 

Table  4 lists  the  ZFA  and  MCD  peaks  for  these  broad  bands 

associated  with  the  y and  y'  systems.  The  vibrational 

spacings  are  close  to  the  values  obtained  for  the  y and  y ' 

transitions.  Due  to  the  breadth  of  the  transitions  the  peak 

maxima  are  uncertain  by  about  20  cm  ^ and  there  is  a shift  to 

lower  energies  in  the  MCD  spectra. 

74 

Recently,  Hildenbrand  has  measured  the  equilibrium 
constant  for  the  gaseous  reaction  Ti  + Ti02  t 2TiO  between 
1975  and  2222°K.  The  equilibrium  constant  was  found  to  vary 
between  5.05  x 10“^  0 1975°K  to  2.72  x 10“^  0 2222°K.  The 
mole  percent  of  TiO  and  Ti  in  a vapor  at  1975”K  is  45%  each 
W'ith  10%  Ti02.  The  percentage  of  Ti02  decreases  slightly  to 

8%  upon  increasing  the  temperature  to  2222°K.  Although 

63  “1 

Weltner  did  observe  infrared  bands  at  940  cm  attributable 

to  Ti02  in  Ne  matrices,  the  effect  of  temperature  on  the 

equilibrium  and  the  opposite  effect  noted  in  observation  5 

above  would  disfavor  Ti02 • By  illumination  with  a mercury 

lamp,  he  also  observed  emission  from  Ti02  in  Ne,  but  the 

O 

( 0 , 0 , 0) - ( 0 , 0 , 0 ) band  of  the  transition  occurred  at  5295  A, 

a much  shorter  wavelength  than  the  present  broad  bands.  The 

fact  that  bulk  Ti02  is  white  would  also  tend  to  eliminate 

7 5 

it  from  consideration.  Electric  deflection  studies  of 
Ti02  also  indicate  that  it  has  a permanent  dipole  moment 
indicative  of  a bent  structure.  Therefore,  the  orbital 
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TABLE  4 


ZFA  and  MCD  Peak  Positions 
Associated  with  the  y and 
TiO  Isolated  in  Argon 


for  Broad  Bands 
Y ' Systems  of 
Matrices 


ZFA 

MCD 

MCD 

A (A) 

A (A) 

V (cm  ) 

7537 

7574 

13203 

Y Band 

7112 

7144 

13998 

6711 

6729 

14861 

6469 

6484 

15423 

Y ' Band 

6134 

6149 

16263 

5834 

5854 

17082 

AG' 


v + 


(cm 


795 

863 


840 

819 
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contribution  to  the  ground  state  magnetic  moment  will  be 
zero  since  it  will  be  orbitally  non-degenerate;  however, 
this  does  not  exclude  a spin  contribution  which  could  yield 
a Cq  term. 

The  close  association  of  the  broad  peaks  with  the  y 
and  y'  systems  indicates  a similarity  between  the  unknown 
species  and  TiO.  The  likely  possibilities  include  aggregates 
of  TiO  molecules  in  a geometry  suitable  for  a ground  state 
magnetic  moment  to  be  generated  or  for  a weak  complex  of  Ti 
and  TiO.  With  the  available  data  it  is  impossible  at  this 
time  to  make  a conclusive  assignment. 

The  y bands  correspond  to  the  transitions  A $-<-X  A^. 

From  the  Boltzman  population  considerations  above  only 

transitions  from  the  ^A^^  level  are  observed  in  these 

studies.  Figure  15  shows  the  ZFA  absorbance  spectrum  of  the 

y band  system.  The  only  MOD  allowed  transition  is  the 
3 3 

A A^.  The  magnetic  field  splitting  is  illustrated  in 

Figure  19,  and  evaluation  of  the  term  which  is  observed 
in  the  MCD  spectrum  is  illustrated  in  Figure  18.  The  site 
splittings  observed  in  the  y ' system  is  much  less  pronounced 
in  the  y system.  Several  matrix  MCD  spectra  at  14 °K  show 
splitting  of  the  negative-going  lobes  of  several  wavenumbers, 
but  these  features  disappear  on  warming  the  matrix  to  24°K 
due  to  the  general  broadening  of  the  transitions  at  higher 
temperatures.  Table  5 lists  the  wavelengths  for  the 
positive-  and  negative-going  lobes  of  the  MCD  of  the  y 
system.  The  zero  crossing  positions  in  Table  5 and  the  ZFA 
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TABLE  5 


MCD  Peak  Positions  for  y System  of 
TiO  Isolated  in  an  Ar  Matrix 


Yl. 

Y (A) 

- / -In 

V (cm  ) 

MCD 

Sign 

0 

7270 

13755 

- 

0 

7235 

13822 

+ 

1 

6853 

14592 

— 

1 

6826 

14650 

+ 

2 

6481 

15430 

- 

2 

6457 

15487 

+ 

3 

6150 

16260 

— 

3 

6133 

16305 

+ 

Zero  (A) 
Crossing 


7253 


6840 


6469 


6142 
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maxima  in  Table  2 agree  within  one  Angstrom  as  would  be 

predicted  for  an  MCD  term. 

The  very  weak  e system  was  first  observed  in  matrices 

by  Weltner  et  al.  and  assigned  to  the  transition  . 

R 

The  only  MCD  allowed  transition  is  illustrated  in  Figure  19 
as  in  the  state.  These  weak  bands  occur  in  the  near  IR 

region  where  instrumental  sensitivity  is  poor  due  to  low 
photomultiplier  tube  response  and  decreasing  light  levels. 
Often,  the  most  sensitive  lock-in  amplifier  range  was 
required  to  detect  the  MCD  signal.  A strong  lamp  emission 
band  occurs  near  the  (0,0)  transition  and  causes  considerable 
difficulty  in  ZFA  spectra  of  that  region. 

The  MCD  spectrum  of  the  (0,0)  band  of  the  e system  is 
shown  in  Figure  22,  and  the  predicted  positive  term  is 
present.  Strong  site  splitting  is  observed  as  in  the  B^n 
state.  Three  distinct  sites  are  noted  in  this  matrix  spec- 
trum, but  in  other  matrices  the  intermediate  energy  site  is 
much  more  dominant. 

The  ultraviolet  region  of  the  spectrum  is  somewhat 

congested  due  to  the  presence  of  many  electronic  transitions 

7 

of  both  TiO  and  Ti.  Pathak  and  Palmer  observed  a triplet- 
triplet  emission  system  designated  D-X^A  in  the  flame  spec- 
trum of  TiO.  They  were  not  able  to  assign  vibrational 
progressions  because  the  spacing  was  irregular.  Shortly 

ft  ^ 

after  this  work  was  published,  Weltner  was  able  to  observe 
this  system  in  absorption  of  TiO  isolated  in  Ne  matrices  at 
4°K.  The  same  type  of  features  were  observed  in  the  matrix 
spectra  which  indicated  that  the  bands  arose  from  at  least 
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two  perturbing  electronic  states  designated  as  systems  I and 
II.  Both  systems  appeared  to  be  perturbed  at  several  vibra- 
tional levels  since  the  vibrational  progressions  do  not 
exhibit  any  regular  behavior. 

An  absorption  spectrum  of  an  argon  matrix  containing 
TiO  and  Ti  is  presented  in  Figure  23.  The  bands  for  Ti 
atoms  are  more  intense  than  TiO  in  this  spectrum  and  corre- 
late well  with  Ar  matrix  spectra  of  Ti  analyzed  by  Gruen  and 

77  . 3 

Carstens.  Since  the  ground  state  of  Ti  is  and  magnet- 

ically degenerate,  strong  terms  are  observed  in  the  MCD 
spectrum  illustrated  in  Figure  24.  These  terms  dominate 
the  weaker  terms  for  TiO,  but  the  temperature  independent, 
derivative  character  of  the  terms  facilitates  the  assign- 
ment of  the  TiO  peaks.  Table  6 presents  the  correlation 

r o 

between  the  Ne  matrix  spectra  of  TiO  obtained  by  Weltner 
and  the  Ar  matrix  spectra  obtained  in  this  work.  The 
presence  of  positive  A^  terms  indicate  that  the  D state  of 

TiO  cannot  be  a state,  but  does  not  distinguish  between 

3 ' 3 7 R 

n and  $ states,  suggested  by  Wentink  and  Spindler  for 

the  D state  of  TiO. 

The  sign  of  the  Cq  terms  observed  for  the  Ti  atomic 

transitions  can  be  used  to  check  the  tentative  assignments 

7 7 

of  Gruen  and  Carstens.  Table  7 lists  the  Ti  atomic  transi- 
tions observed  in  Ar  matrices  along  with  the  observed  MCD 
signs  and  calculated  MCD  Cq  terms  based  upon  the  proposed 
assignments.  Two  of  the  MCD  signs  are  in  conflict  with  the 
assignments  of  Gruen  and  Carstens.  The  band  at  2874  A has 
a negative  MCD  Cq  term,  but  the  value  predicted  from  their 
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TABLE  6 


Ultraviolet  Transitions  Observed 
for  TiO  in  Argon  and  Neon  Matrices 


Ne 


Matrix* 

^ (X) 

Ar  Matrix^ 
A (A) 

MOD  Shape 

System 

3117 

3137 

+A^  term 

I 

3061 

3087 

+A^  term 

II 

3013 

3032 

+A^  term 

I 

2969 

2998 

+A^  term 

II 

2932 

2956 

+A^  term 

I 

2879 

2907 

+A^  term 

II 

2805 

2831 

+A^  term 

II 

* Reference  63. 
t This  study. 


Ti  Atomic  Transitions  Observed  in  an  Argon  Matrix  at  14 °K 
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3 . ... 

assignment  is  positive.  The  negative  MCD  term  indi- 
cates that  this  transition  is  to  the  state  observed  in 

gas  phase  spectra  at  2934  A.  The  strong  peak  at  2550  A 
displays  a strongly  positive  term  which  would  be  forbidden 
if  the  ^F^  assignment  were  correct.  The  MCD  indicates  that 
the  state  is  probably  the  state  located  at  2633  A in 

the  gas  phase  spectra. 

Three  other  medium  intensity  bands  are  observed  at 
2985  A,  2842  A and  2808  A,  The  two  longer  wavelength  bands 
have  positive  terms,  and  the  other  band  has  a negative 
Cq  term.  These  transitions  were  not  observed  by  Gruen  and 
Carstens  but  are  probably  due  to  Ti  atoms. 


CHAPTER  VI 


RESULTS  OF  ZrO  AND  HfO  EXPERIMENTS 
Introduction 

Due  to  the  effects  of  the  lanthanide  contraction,  both 

the  atomic  radii  of  Zr  (1.45  A)  and  Hf  (1.44  A)  and  the 

radii  of  the  ions  Zr’’"'*  (0.74  A)  and  Hf'^'^  (0.75  A)  are  nearly 

identical.  The  chemical  properties  of  the  compounds  formed 

from  these  elements  are  more  similar  than  any  other  pair  of 

cogeneric  elements.  In  all  geological  deposits  of  Zr,  Hf  is 

found  in  the  range  of  fractions  of  a percent,  and  their 

48 

separation  is  very  difficult.  As  mentioned  previously, 

the  ground  states  of  ZrO  and  HfO  are  the  same,  and 

different  from  the  ground  state  of  the  other  group  IVb  oxide 
3 

TiO,  A . The  reversal  of  ground  state  configuration  of 

K 

2 

(a6)  for  TiO,  as  shown  in  Figure  14,  to  (a)  for  ZrO  and  HfO 

is  presumably  due  to  a combination  of  heavy  atom  properties 

37 

and  electron  correlation  effects.  This  chapter  involves  a 
discussion  of  the  experimental  results  obtained  for  ZrO  and 
HfO  isolated  in  Ar  matrices. 

Results  for  Zirconium  Monoxide 

Like  TiO,  ZrO  has  been  studied  for  years  because  of  its 
astrophysical  importance  in  Class  S stars,  e.g.,  R And  and 
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R Cyg.  One  of  the  characteristic  differences  between  the 
low  surface  temperature  R and  S stars  is  that  in  S stars  the 
bands  for  ZrO  are  more  prominent  than  TiO  bands.  Two  ZrO 
singlet  bands,  designated  A and  B,  were  first  identified  in 

on  o 1 

gas  phase  spectra  taken  in  1950,  ' and  six  triplet 

emission  systems  have  also  been  observed.  In  1957, 

o .82 

Akerlind  showed  that  the  A and  B transitions  do  not  involve 
a common  upper  or  lower  state.  The  first  matrix  isolation 

ft  ') 

study  was  carried  out  by  Weltner  and  McLeod  in  1965. 

Until  their  work,  the  ground  state  of  ZrO  was  considered  to 
3 

be  A like  TiO,  but  since  they  only  observed  transitions 
originating  in  the  state  at  very  low  temperatures,  the 

ground  state  designation  was  changed  to 

The  ZrO  matrix  deposition  conditions  used  in  this  study 
were  those  used  by  Weltner  and  McLeod.  The  Knudsen  cell  was 
loaded  with  Zr02  (99.9%  Alfa),  and  deposition  carried  out  at 
2630°K. 

Table  8 lists  the  wavelength  absorption  bands  of  the 

r o 

ZFA  observed  for  ZrO  in  Ar  matrices,  by  Weltner  and  McLeod 
along  with  the  MOD  and  ZFA  peak  positions  found  in  this 

8 3 

study.  Several  ZrO  bands  were  observed  by  Bates  and  Gruen 
in  a recent  study  of  Zr  atoms  isolated  in  rare  gas  matrices, 
and  these  are  also  listed  in  Table  8.  It  should  be  noted 
that  the  wavelength  of  the  peak  observed  by  Bates  and  Gruen 
at  5681  A and  attributed  by  them  to  the  (3,0)  band  of  ZrO  is 

. ® O 

shifted  by  22  A from  the  peak  at  5659  A.  A state  diagram 
for  the  known  transitions  and  states  of  ZrO  is  presented  in 
Figure  25. 
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TABLE  8 


ZFA 

and  MCD 

Bands 

for  ZrO 

Gas  Band 

ZFA* * 

/ 

System 

v' 

A (A) 

ZFA 

ZFA*^ 

MCD*^ 

Sign 

C 

0 

6633 

6634 

6635 

6674 

6598 

+ 

1 

6270 

6276 

6275 

6301 

- 

6244 

+ 

2 

5947 

5951 

5950 

5971 

- 

5930 

+ 

3 

5659 

5681 

5659 

5675 

- 

5647 

+ 

4 

5406 

5400 

Too  Weak 

D 

0 

6024 

6021 

6024 

6020 

+ 

1 

5729 

5728 

5728 

+ 

2 

5461 

5471 

5460 

Too  Weak 

Forbidden 

0 

5302 

5304 

5320 

- 

Transition 

5292 

+ 

1 

5084 

5080 

5098 

- 

5064 

+ 

A 

0 

3678 

3678 

1 

3572 

3570 

No  MCD 

2 

3470 

3469 

3371 

3260 

^ o 

* Results  of  Weltner  and  McLeod. 

8 3 

t Results  of  Bates  and  Gruen. 

/ Results  of  this  study. 
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The  band  system  with  its  (0,0)  transition  at  6635  A 
in  Ar  is  designated  as  the  C system  and  corresponds  to  the 
gas  phase  singlet  system  at  6495  A.  Weltner's  assignment 

of  the  transition  to  e^n-<-X^E^  was  confirmed  by  Balfour  and 

84  . 

Tatum  in  an  analysis  of  the  arc  emission  spectrum  of  ZrO. 
This  region  of  the  emission  spectrum  is  dominated  by  the 
Y system,  A^<J>-X^A.  The  red-degraded  head  of  the  C system 
occurs  between  the  (0,0)  and  (1,1)  band  heads  of  the  y 
system  and  is  completely  overlapped  by  the  y system  which 
makes  the  interpretation  of  this  region  very  difficult. 

The  predicted  magnetic  field  splittings  and  circularly 
polarized  transitions  for  the  C system  of  ZrO  are  illustrated 
in  Figure  26.  Since  the  ground  state,  is  magnetically 

non-degenerate,  the  ground  state  will  not  be  split,  and  the 
only  possible  MCD  terms  for  ZrO  are  and  8^  terms.  The 
predicted  MCD  term  for  the  C system  is  +13  and  is  readily 

confirmed  by  the  positive  temperature  independent  A^  term 
observed  for  this  band  in  Figure  27.  The  ZFA  spectrum  for 

this  wavelength  region  is  also  shown  in  Figure  27. 

6 2 

Weltner  noted  weak  matrix  site  effects  in  the  wings 
of  the  ZrO  bands  in  Ar.  Small  amounts  of  blue  shading 
occurs  in  the  ZFA  bands  of  the  C system.  Weak  site  effects 
are  also  noted  to  the  blue  of  the  (0,0)  and  (1,0)  and  as  a 
broadening  of  the  (2,0)  and  (3,0)  bands  in  the  MCD.  The 
intensity  of  these  sites  changed  on  annealing  the  matrix. 

The  spacing  between  the  sites  on  the  (0,0)  and  (1,0)  MCD 
bands  is  851  cm  ^ which  is  very  close  to  the  spacing  of 
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Figure  26. 


Magnetic  field  splitting  and  MCD  allowed 
transitions  of  ZrO  and  HfO. 


Q 

'V 

C • 

(I3  X 
•H 
U M 
-p 
(U  fC 

^ e 

-p 

u 

iw  <C 

o 

C 

(C 

p 

-P  c 
U -P 

<D 

CUO 
tn  p 
N] 

C 

fa  tp 
N O 


T3  W 

c g 

(0  (U 
■P 

Q cn 
U >1 
S tfl 


(N 

0 

P 

tn 

•rH 

fa 


I 


aoupqjosqv  V 


115 


859  cm  ^ between  the  nearby  positive  lobes  of  the  (0,0)  and 
(1,0)  MCD  bands. 

The  (2,0)  band  of  the  C system  overlaps  with  the  (0,0) 
band  of  the  D system.  In  the  ZFA  spectrum  shown  in  Figure  27 
the  intensity  of  the  bands  for  the  D system  is  greater  than 
the  overlapping  bands  belonging  to  the  C system  at  shorter 
wavelengths,  but  in  the  MCD  spectrum  of  Figure  27  this 
situation  is  reversed  and  the  C system  bands  are  relatively 
more  intense.  Both  MCD  band  systems  show  no  variation  with 
temperature,  and  the  bands  for  the  D system  do  not  show  the 
characteristic  derivative  shape  for  an  term  but  follow 
the  lineshape  function.  The  maximum  of  the  D system  bands 
in  the  ZFA  and  MCD  spectra  correspond. 

The  assignment  of  the  D system  to  the 

^n(7T^6)  X^E^(Tr'^a^)  transition  was  deemed  very  uncertain  by 

6 2 

Weltner  and  could  have  also  been  assigned  to  the 
^Z^(ao*)  X^Z^(tt‘^o^)  transition.  Both  of  these  transitions 

would  be  allowed  in  ZFA  spectra.  The  transition  to  the 
state  should  yield  a positive  term  in  MCD  spectra; 
whereas  the  transition  to  the  state  would  be  forbidden 

in  MCD  unless  there  was  a Bq  term  induced  by  magnetic  field 
mixing  with  a nearby  state.  The  most  likely  candidate  for 
this  mixed  in  state  is  the  C system  in  which  the  e^II  state 
does  possess  a magnetic  moment.  The  assignment  of  the 
D system  to  the  ^Z^(oa*)  X^Z^(Tr^a^)  transition  seems  more 

logical  in  light  of  the  MCD  sign,  shape,  intensity,  tempera- 
ture independence,  and  the  presence  of  the  nearby  e^H  state 

8 5 

for  mixing  to  form  a Bq  term.  Recently  Phillips  and  Davis 
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have  examined  this  wavelength  region  closely  in  the  arc 
emission  spectrum  of  ZrO  and  discovered  a 

transition  with  the  (0,0)  band  at  5860  K,  which  is  close  to 
the  6024  A band  in  Ar  matrices.  The  spacing  between  the 

vibrational  levels  is  870,  867  and  855  cm  ^ for  ^^ZrO,  ^^ZrO 

96  — 1 

and  ZrO,  respectively,  which  is  also  close  to  the  858  cm 

spacing  for  the  vibrational  levels  of  the  D system  in  Ar 

matrices . 

An  example  of  a ^Z^  ^Z^  transition  where  there  is  no 

nearby  state  available  to  form  a 8q  term  in  MCD  is  demon- 
strated in  the  A system  of  ZrO.  This  system  shows  a 
prominent  progression  in  ZFA  spectra  in  the  range  between 
3700  A and  3200  A as  illustrated  in  Figure  28.  Weltner 
noted  three  members  of  the  progression,  and  two  more  have 
been  added  in  this  study.  The  spacings  between  the  vibra- 
tional bands  are  not  quite  regularly  anharmonic,  but  this 
may  be  due  to  overlap  with  Zr  atomic  bands  which  overlap 
the  3371  A and  3678  A bands.  The  Zr  atomic  transitions  are 
discussed  below. 

The  original  assignment  of  the  A system  by  Uhler  and 

81  11 
Akerlind  was  doubtful  because  a Z-  Z transition  could 

not  be  distinguished  from  a transition.  Weltner 

favored  the  former  assignment  on  the  basis  of  the  corre- 
spondence between  the  infrared  vibrational  frequency  of 
ZrO  in  a Ne  matrix  (975  cm  and  that  in  the  lower  state 
of  the  A gas  emission  system  (969.7  cm  , rather  than  in 
the  lowest  ^A  state  (931  cm  . The  results  of  this  MCD 
study  further  confirm  this  assignment.  The  MCD  spectrum 
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of  the  wavelength  region  of  the  A system  (see  Figure  29) 
shows  only  those  transitions  attributable  to  Zr  atomic  Cq 
terms  and  none  for  the  MCD  forbidden  transition  of 

the  A system. 

O 

Two  very  weak  transitions  are  observed  at  5304  A 

6 2 

and  5080  A which  were  also  observed  by  Weltner  and 

3 1 + 

assigned  by  him  to  the  forbidden  transition  B II  X Z , 
based  on  an  upper  state  vibrational  frequency  of  836  cm  ^ 
in  Ne  matrices  as  compared  with  the  known  gas  phase 
v'  = 838.8  cm  ^ for  the  state.  Although  the  AS=0 

selection  rule  is  violated  in  transitions  to  all  three  of 

3 ■ 3 

the  n multiplet  states,  only  the  transition  to  the  B 

state  does  not  also  violate  the  AZ=0  selection  rule  as 

illustrated  in  Figure  30.  Note  also  that  the  MCD  transition 

to  the  JI=1  multiplet  is  the  only  one  that  does  not  violate 

the  AM^=0  selection  rule  for  MCD.  Lauchlan,  Brom  and 

ft  6 . ° . + 

Broida  excited  this  state  using  the  5145  A line  of  an  Ar 

ion  laser  on  a Ne  matrix  containing  ZrO  where  the  transition 

occurred  at  5154  A.  They  observed  strong  B^II^  X'^A^ 

emission  with  this  excitation  which  strongly  supports  the 

assignment.  Weak  MCD  transitions  are  observed  for  this 

system  (see  Figure  31)  with  positive  terms  which 

further  supports  the  assignment  since  the  predicted  term 

for  the  B^II^  X^Z^  is  also  positive. 

6 2 

Five  other  weak  bands  were  noted  by  Weltner  for  ZrO 
matrices  in  Ar,  and  four  of  these  bands  can  now  be  assigned. 
Two  of  these  weak  bands  at  4855  A and  4632  A show  positive 
A^  terms  in  the  MCD  (see  Figure  31) . Their  vibrational 
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Figure  30.  Magnetic  field  splittings  and 
terms  for  the^n-<-^E''‘  transition 
of  ZrO. 
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spacing  is  992  cm  with  some  site  splitting  in  the  MCD 
shown  in  Figure  31.  Unfortunately,  only  spectra  of  Zr^^O 
were  obtained  in  Ar  matrices,  and  therefore,  any  shifts 
that  might  help  assign  the  (0,0)  transition  for  these 
bands  are  not  known.  Zr  0 spectra  were  obtained  by  Weltner 
for  Ne  matrices;  however,  these  transitions  were  not 
noted  in  the  Ne  matrix  spectra.  This  transition  may  be 
the  ^n(7T^6)  transition  previously  tentatively 

assigned  to  the  D system  which  was  revised  to  the 
L (ao  ) ^ X Z (tt  a ) , vvde  supra. 

Two  of  the  other  weak  transitions  observed  in  ZFA 

spectra  in  Ar  matrices  at  3770  A and  3721  A are  readily 

attributable  to  Zr  atomic  transitions  which  were  observed 

8 3 

recently  by  Bates  and  Gruen.  These  transitions  are 

readily  observed  in  the  MCD  spectra  since  their  overlap 

with  the  A bands  of  ZrO  is  removed  by  the  forbidden  character 

of  the  MCD  of  the  A bands  (see  Figure  29) . The  ground  state 
3 

of  Zr  xs  F^  and  therefore  shows  strong  term  behavior 
in  MCD  spectra.  Although  the  experiments  in  this  study 
were  directed  toward  ZrO,  enough  Zr  atoms  are  present  under 
these  experimental  conditions  to  observe  all  of  the  moder- 
ately intense  Zr  atomic  transitions  observed  by  Bates  and 
Gruen.  These  transitions  are  listed  in  Table  9.  Although 
Bates  and  Gruen  did  not  make  any  correlations  between  the 
gas  phase  spectra  and  the  observed  matrix  transitions,  the 
sign  of  the  MCD  terms  for  the  Zr  atomic  transitions  allows 
prediction  of  the  values  for  the  upper  states  which  are 
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TABLE  9 

Zr  Atomic  Transitions  Observed  in  Ar  Matrices 


0 + 
A (A)^ 

Strength 

o J 

A (A)*^ 
This  Work 

MCD 

Sign 

Final 

9. 

* 

4494 

3 

4367 

m 

4368 

- 

3 

4300 

vw 

4322 

+ 

1 or  2 

4005 

vw 

4028 

+ 

1 or  2 

3819 

m 

3825 

- 

3 

3771 

vs 

3770 

- 

3 

3719 

s 

3721 

+ 

1 or  2 

3459 

w 

3457 

- 

3 

3401 

vs 

3400 

- 

3 

3349 

s 

3350 

+ 

1 or  2 

3050 

m 

3037 

+ 

1 or  2 

2906 

m 

2901 

+ 

1 or  2 

2865 

s 

2848 

+ 

1 or  2 

2812 

vw 

2796 

- 

3 

2729 

w 

2733 

+ 

1 or  2 

2691 

ms 

2689 

+ 

1 or  2 

2594 

w 

2600 

3 

* Not  reported. 

n o 

t Results  of  Bates  and  Gruen. 

/ Results  of  this  study. 
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also  presented  in  Table  9.  All  of  the  peaks  listed  show  the 

expected  temperature  variations  for  an  MCD  Cq  term. 

o ° . 

The  ZFA  spectrum  from  2500  A to  3250  A is  shown  in 

Figure  32.  A large  number  of  strongly  absorbing  bands  are 

noted  in  this  region,  some  of  which  are  attributable  to  the 

Zr  atomic  transitions  listed  in  Table  9.  The  majority  of 

the  bands  in  this  region  are  not  due  to  Zr  atoms,  but  are 

probably  due  to  previously  unidentified  bands  of  ZrO.  The 

6 2 

shortest  wavelength  bands  noted  by  Weltner  comprise  the 

O ^ . . 

A system  between  3700  A and  3200  A.  Although  no  mention  is 

made  for  the  shortest  wavelength  examined  in  the  ZrO  spectra, 

he  does  note  that  the  region  below  3100  A was  not  examined 

8 0 

in  the  HfO  spectra.  In  the  arc  emission  studies  of  Afaf 
in  1950,  the  shortest  wavelength  bands  observed  was  the  (p 
triplet  system  between  3121  A and  2940  A with  the  (0,0) 
transition  located  at  2941  A.  This  system  was  very  weak, 
and  only  transitions  from  the  v'=0  level  were  observed. 

The  low  intensity  of  emission  is  not  surprising  since  this 
state  is  over  34000  cm  ^ above  the  ground  state  and  the 
fractional  population  in  the  excited  state  is  very  low. 

The  MCD  spectrum  of  this  region  is  shown  in  Figure  33. 
The  Zr  atomic  bands  are  much  enhanced  in  the  MCD  spectra 
since  they  are  due  to  Cq  terms  which  are  much  more  intense 
than  terms.  An  example  of  this  enhancement  can  be  seen 

O 

by  comparing  the  small  shoulder  Zr  transitions  at  2689  A 
in  the  ZFA  spectrum  of  Figure  32  with  the  much  stronger  peak 
in  the  MCD  spectrum  of  Figure  33. 
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Figure  33.  MCD  spectrum  of  the  ultraviolet  region  of  an 
Ar  matrix  containing  Zr  and  ZrO. 


132 


Table  10  lists  the  peak  positions  for  the  ultraviolet 

ZrO  bands  in  the  ZFA  and  MCD  spectra.  In  all  cases  the 

middle  of  the  ZFA  band  is  located  in  the  middle  of  the 

two  MCD  lobes.  This  agreement  is  readily  seen  in  Figure  34 

where  both  spectra  are  plotted  with  the  same  abcissa.  Each 

MCD  band  appears  as  a temperature  independent,  positive 

term  which  indicates  that  the  transitions  are  most 

likely  More  than  one  transition  is  indicated  by 

the  Deslandres  table  illustrated  in  Figure  35.  States  that 

have  overlapping  energy  levels  often  perturb  each  other  so 

that  the  vibrational  progressions  do  not  show  a regular 

anharmonic  behavior  as  seen  in  the  TiO  ultraviolet  bands  of 

System  I and  II.  A possible  grouping  of  vibrational  bands 

is  indicated  in  Table  10;  however,  due  to  the  strong 

overlap,  these  assignments  are  very  tentative. 

The  possibility  that  these  bands  are  due  to  Zr02  is 

not  completely  eliminated.  Weltner  observed  one  strong 

band  at  965  cm  ^ for  Zr^^O  in  Ar  matrices  and  other  bands 

— 1 8 7 

at  887,  880,  822,  816  and  794  cm  . Linevsky  observed 

bands  at  819  and  884  cm  ^ in  the  infrared  spectra  of  ZrO 

in  Ar  matrices  which  he  attributed  to  ZrO^  absorption.  The 

observed  frequencies  differ  by  about  5 cm  ^ between  Weltner 's 

work  and  Linevsky ' s data.  In  mass  spectrometer  studies  of 

8 8 

the  vaporization  of  Zr02,  Chupka,  Berkowitz  and  Inghram 

observed  Zr02  as  a component  in  the  vapor  state,  although 

ZrO  was  the  dominant  species  with  the  ( ZrO) ^/ ( Zr02 ) ^ ratio 

89 

of  about  6.  Conlon  and  Doyle  have  examined  the  optical 
absorption  of  thin  films  of  Zr02  in  the  region  between 


33389 

33956 

34364 

34771 

35273 

35486 

35881 

36232 

36900 

37793 

38314 

38910 

39557 
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TABLE  10 


Ultraviolet  Bands  Observed  for  ZrO  in  Ar 


Temperature 


Abs 

MCD 

Sign 

Dependence 

2995 

2997 

None  - 

2993 

+ 

2945 

2951 

— 

None 

975  cm  ^ 

2943 

+ 

2910 

2916 

— 

None  - 

2895 

+ 

2876 

2880 

— 

None 

909  cm  ^ 

2867 

+ 

2835 

2840 

Nnnf^  —I 

2830 

+ 

2818 

2822 

- 

None 

959 

cm 

2812 

+ 

2787 

2790 

— 

None  1 

2784 

+ 

2760 

2765 

— 

None 

1019  J 

cm 

2754 

+ 
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2050  and  10000  A.  These  films  were  generated  by  oxidation 
of  thin  films  of  Zr  deposited  on  fused  quartz  disks. 

Therefore  their  spectra  are  properties  of  solid  state 
The  absorption  threshold  of  Zr02  is  about  2300  A and  rises 
sharply  to  their  limit  of  measurement  at  2050  A.  These 
authors  attribute  the  absorption  to  exciton  formation  or 
band-to-band  transitions  which  would  probably  not  be  observed 
in  Ar  matrices. 

If  the  transitions  were  due  to  Zr02,  the  only  way  the 
MCD  would  be  seen  is  if  Zr02  were  linear  so  that  orbital 
degeneracy  could  occur  or  if  the  transitions  were  forbidden 
assuming  Zr02  has  a singlet  ground  state.  These  possibili- 
ties do  not  seem  likely  since  Ti02  has  been  shown  to  be 
bent,  and  the  intensity  of  the  transitions  seems  to  rule 
out  intermultiplicity  transitions. 

The  group  of  bands  in  the  region  3025  A to  3140  A is 
observed  in  spectra  of  ZrO  and  HfO  (see  Figure  36) . These 
bands  are  much  less  intense  in  the  ZFA  spectra  (see 
Figure  32)  and  are  readily  attributable  to  the  R system  of 

TaO  as  observed  in  the  gas  phase  at  3088  A by  Cheetham  and 
90 

Barrow.  These  bands  were  also  observed  in  our  studies  of 

3 8 

TaO  isolated  in  Ar  matrices.  The  other  bands  observed  in 
matrices  of  ZrO  and  less  intensely  in  HfO  matrices  are 
listed  below: 
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The  presence  of  TaO  in  the  matrices  is  due  to  oxidation  of 
the  Ta  end-caps  in  the  Knudsen  cell  at  the  high  temperatures 
required  for  formation  of  HfO  and  ZrO. 

Results  of  Hafnium  Monoxide 

A scant  five  years  after  the  discovery  of  the  element 
91 

hafnium.  Meggers  published  the  first  study  of  HfO  in  1928. 

92 

This  report  was  followed  by  the  work  of  King  in  1929  and 

9 3 

Shaw  and  Ketcham  in  1934;  however,  no  detailed  level 

assignments  could  be  made  at  that  time.  In  1951, 

94 

Krishnamurti  published  a preliminary  report  on  HfO  but 

never  followed  up  with  a detailed  analysis.  The  Vatican 
95 

Atlas,  published  in  1957,  gave  spectra  and  preliminary 
assignments  for  9 sequences  in  HfO. 

The  first  matrix  isolation  study  of  HfO  was  published 
in  1965  by  Weltner  and  McLeod. By  using  ^^O  labeling 
techniques,  they  were  able  to  determine  the  (0,0)  transitions 
for  six  band  systems  in  Ar  and  Ne  matrices  at  4 and  20°K. 

The  wavelengths  obtained  for  the  bands  of  these  systems  are 
presented  in  Table  11  along  with  ZFA  and  MOD  spectral  wave- 
lengths obtained  in  the  present  study  under  the  same  matrix 
deposition  conditions. 

As  noted  above,  the  properties  of  Zr  and  Hf  are  very 
close,  and  this  similarity  extends  to  the  oxides.  The 
ground  state  constants  for  ZrO  and  HfO  are  listed  on  the 
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Term 

Reduced 

Mass 

Tp  (A) 

Wp  (cm  ) 

ZrO 

(X^Z) 

13.55 

1.711 

969.8 

<1 

X 

1.728 

936.5 

HfO 

(X^Z) 

14.69 

1.723 

974.1 

The  correlation  of  co^  values  in  the  ground  state  carries 

over  into  the  excited  states  where  the  ZrO  bands  shift  very- 

little  in  going  to  HfO  bands.  This  correspondence  holds 

very  well  in  Ne  matrices  but  less  well  in  Ar  matrices,  as 

shown  in  Table  11,  and  was  used  by  Weltner  and  McLeod  to 

assign  the  D,  E and  F states  as  ^II-^^Z,  and  ^Z-<-^Z, 

respectively.  These  assignments  agree  with  the  designations  by 

95 

Gatterer  et  al.  based  on  the  appearance  of  the  gas  phase 

bands.  The  A,  F and  G systems  have  the  appearance  of  ^Z^^Z 

transitions,  while  the  B,  D and  E systems  are  typical  of 

^n-^Z  transitions.  This  gas  phase  analysis  was  further 

96 

refined  by  Edvinsson  and  Nylen  in  1971  in  an  examination 

of  13200  spectral  lines  of  the  A,  B,  D,  E,  F and  G systems. 

97 

Later,  Edvinsson  revised  the  assignment  of  the  A and  B 

1 3 1 3 

systems  from  A Z to  the  and  B IT  to  the  levels  of 

a triplet  state.  The  high  resolution  data  fit  either 

assignment. 

The  ZFA  and  MCD  spectra  of  the  two  ^Z-^Z  ultraviolet 
bands  F and  G are  presented  in  Figure  37.  The  ZFA  shows  a 
distinct  3 member  progression  for  both  systems,  but  the  MCD 
shows  no  signal  as  expected  for  transitions  between  non- 
degenerate states. 

. O O 

The  spectral  region  between  3700  A and  4400  A contains 
the  overlapping  D and  E systems.  The  ZFA  spectrum  of  the 


143 


TABLE  11 

Correlation  between  ZrO  and  HfO  Bands 
in  Ne  and  Ar  Matrices 


Compound 

Band 

System 

Ne 

Matrix 
A (A) 

AG  , , 
(cm  ^ ) 

Ar 

Matrix 

A (A) 

AG  , , 
(cm  A ) 

ZrO 

A 
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838 

3678 

807 
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F 

3642.5 

838 
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833 
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854 
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region  is  shown  in  Figure  38,  and  the  MCD  spectrum  is  in 
Figure  39.  The  gas  phase  analysis  of  Edvinsson  and  Nylen^^ 
shows  that  there  are  perturbations  between  the  D(v'=2)  and 
E(v'=0)  levels,  the  D(v'=3)  and  E(v'=l)  levels,  etc.  These 
perturbations  occur  in  the  regions  where  the  stronger 
D system  overlaps  the  E system.  The  perturbations  are 
manifest  by  irregular  isotopic  effects  in  the  0-0  sequence 
of  the  E system  and  by  the  non-proportionality  between 
J(J+1)  and  the  A-type  doubling  in  perturbed  levels.  For 
the  E^n(v'=l)  level,  the  A-type  doubling  changes  sign  at 
J=106. 

Both  D and  E systems  are  assigned  to  states  which 
show  magnetic  field  splittings  as  illustrated  in  Figure  26. 
The  predicted  positive  term  is  readily  evident  in  the 
D system  bands.  The  less  intense  bands  for  the  E system 
are  broad  and  overlap  with  the  D bands,  which  makes  the 
terms  less  evident.  The  middle  of  the  derivative  shape 
for  the  E bands  in  MCD  correspond  to  the  band  maxima 
observed  in  the  ZFA  spectra,  which  indicates  that  the 
assignment  for  the  E state  of  ^11  is  correct. 

The  A and  B band  systems  of  HfO  do  not  have  any 

counterparts  in  the  ZrO  spectra;  therefore,  Weltner  and 
6 2 

McLeod  were  unable  to  assign  these  states.  They  presumed 
that  they  are  forbidden  transitions  to  lower  triplet  states 
which  are  made  more  "allowed"  by  the  internal  heavy  atom 

effect  in  the  heavier  HfO  molecule.  The  analysis  of  Gatterer 

95  . 1 

et  al.  assigned  the  A system  to  a Z state  and  the 

1 96 

B system  to  a n state.  Edvinsson  and  Nylen 
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Figure  38.  ZFA  spectrum  of  the  (^n-^Z)  D and  E bands  of 
HfO  in  an  Ar  matrix. 
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MCD  spectrum  of  the  D and  E bands 

of  HfO  in  an  Ar  matrix. 


Figure  39. 
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with  this  assignment  based  on  analysis  of  very  high  resolu- 
tion spectra  obtained  using  a Jarrel  Ash  4.5m  Fastie 
spectrograph  with  a dispersion  between  3.5  and  6.5  A/mm  and 
a microwave  discharge  source.  These  level  assignments  were 
later  changed  by  Edvinsson.  The  high  resolution  data  fit 

either  assignment,  but  Edvinsson  attributes  the  A and  B sys- 
3 3 

terns  to  and  11^  levels  of  a triplet  state,  respectively. 

Intermultiplicity  transitions  are  favored  in  heavier 
molecules  such  as  HfO,  but  the  spin-orbit  coupling  constant 
for  Hf  is  1300  cm  ^ whereas  the  spacing  between  the  band 
heads  for  the  A and  B systems  is  935  cm  This  lower 
energy  separation  would  tend  to  disfavor  the  assignment  as 
multiplets  of  a ^11  state. 

The  MCD  transitions  from  the  ground  state  to  the 
3 3 

and  IIq  states  have  predicted  positive  and  negative 
terms,  respectively.  The  transition  to  the  state  is 

allowed  in  Hund ' s case  a coupling  or  in  case  c coupling 
which  would  be  expected  to  be  favored  by  the  high  molecular 
weight  of  the  molecule.  The  transition  to  the  ^11^  state  would 
not  be  allowed  in  either  case  a or  case  c (see  Figure  40) . 

The  ZFA  and  MCD  spectra  of  the  A and  B bands  are 
presented  in  Figure  41.  The  intensity  of  these  bands  is 
very  weak,  and  they  are  also  rather  broad.  The  MCD  spectrum 
was  obtained  on  the  most  sensitive  scale  of  the  lock-in 
amplifier,  and  the  signal- to-noise  ratio  is  not  as  high  as 
on  other  spectra  obtained  in  this  study.  The  observation 
that  the  A and  B bands  have  different  shapes  in  both  ZFA  and 
MCD  spectra  would  disfavor  the  intermultiplicity  assignment. 


Figure  40.  Predicted  magnetic  field  splitting 
and  terms  for  the  forbidden 

transitions  and  of 

HfO. 
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It  is  interesting  to  note  that  the  A bands  at  5957  A and 
5944  K are  reversed  in  intensity  from  Weltner  and  McLeod 's^^ 
study.  In  their  work  they  observed  a reversible  increase  in 
intensity  of  the  5944  A band  on  cooling  the  matrix  from  20°K 
to  4.2°K  and  argued  that  this  was  due  to  site  effects.  The 
site  effect  hypothesis  is  supported  by  the  reversal  of 
intensity  observed  in  this  study  which  is  obstensibly  due  to 
different  matrix  deposition  conditions.  Weltner  and  McLeod's 
matrices  were  deposited  at  4.2°K,  whereas  14 °K  was  the 
deposition  temperature  in  this  study. 

The  shape  of  the  MCD  of  the  A band  lends  support  to  the 
earlier  assignment  as  since  no  derivative  shape  is 

observed,  which  would  be  expected  for  a The  presence 

of  the  negative-going  MCD  band  is  attributable  to  a Bq  term 
due  to  mixing  with  the  nearby  IT  state  of  the  B system. 

By  use  of  labeling,  Weltner^^  observed  two  (0,0) 

bands  (5589  and  5649  A)  in  the  B system  which  overlap  each 
other.  The  shape  of  the  MCD  bands  do  not  conclusively 
indicate  the  presence  of  a singlet  or  triplet  n state. 

Several  bands  are  observed  in  the  ZFA  spectra  of  the 
ultraviolet  region  between  2500  A and  2900  A.  These  bands 
are  listed  in  Table  12  and  the  spectra  are  presented  in 
Figure  42.  The  four  bands  observed  between  2650  A and 

O 

2900  A are  due  to  TaO  present  in  the  matrix.  The  intensity 
ratio  of  these  peaks  to  the  HfO  bands  is  not  constant  for 
different  matrices,  and  the  wavelengths  match  those  found 


for  TaO  matrices. 
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TABLE  12 

ZFA  Peak  Positions  in  the  Ultraviolet 
Spectra  of  HfO  in  Ar  Matrices 


ZFA 

A (A)  Comment 
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2683 


TaO 

TaO 

TaO 

TaO 


2627 

2613 

2603 


+ term  in  MCD 
+ term  in  MCD 
+ A^  term  in  MCD 
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The  group  of  three  peaks  near  2600  A all  display  +A^ 
term  behavior  in  the  MCD  and  are  probably  due  to  a band 
previously  observed  by  Edvinsson  and  Nylen^^  in  gas  emission 
studies  at  2735.5  A.  This  correlation  implies  a matrix 
shift  to  the  blue  of  about  1700  cm  The  presence  of  the 

positive  term  implies  that  this  state  is  ^n. 

Unfortunately,  in  the  experiments  of  Edvinsson  and 
Nylen,  the  region  between  2300  and  2700  A was  completely 
obscurred  by  SiO  bands,  and  no  correlation  can  be  made 
between  the  matrix  bands  observed  below  2600  A and  gas 
phase  spectra. 


CHAPTER  VII 


CONCLUDING  REMARKS 

The  combination  of  magnetic  circular  dichroism  and 
matrix  isolation  techniques  has  been  shown  by  this  study  to 
be  a powerful  tool  for  the  molecular  spectroscopist . By 
employing  an  asymmetric  probe  such  as  circularly  polarized 
light  and  the  unique  selection  rules  for  MCD,  it  is  possible 
to  determine  more  readily  the  electronic  symmetry  of  molecu- 
lar species  than  in  ZFA  spectra.  The  use  of  the  various 
lineshape  and  temperature  dependent  parameters  of  MCD  gives 
the  spectroscopist  more  criteria  to  use  in  the  assignment  of 
molecular  and  atomic  states. 

The  future  of  this  technique  is  bright  and  should  be 
rather  productive  since  so  few  species  have  been  studied 
thus  far.  As  illustrated  in  the  spectra  of  TiO,  MCD  will  be 
very  helpful  in  studies  of  matrix  site  effects  since  MCD  is 
able  to  resolve  much  finer  details  than  have  been  observed 
previously  in  ZFA  spectra.  The  use  of  MCD  to  determine 
magnetic  moments  is  also  a field  which  should  prove  very 
fruitful . 

Several  possibilities  exist  for  improvement  of  the 
apparatus  built  for  this  experiment.  The  results  have 
shown  that  a relatively  modest  magnetic  field  (5500  gauss) 
is  sufficient  to  obtain  excellent  MCD  data;  however,  the 
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small  diameter  and  long  axial  holes  in  the  pole  pieces  of 
standard  electromagnets  precludes  circularly  polarized 
emission  studies.  A more  compact  design  split  coil,  super- 
conducting magnet  or  a high  field,  rare  earth  permanent 
magnet  should  improve  the  chances  for  observation  of  circu- 
larly polarized  luminescence. 

The  size  of  the  PET  microcomputer  memory  should  be 
expanded  to  32K  in  order  to  be  able  to  store  a full  8000  to 

O 

2500  A scan  in  memory.  This  would  eliminate  the  necessity 
of  storage  of  partial  spectra  on  magnetic  tapes  for  subse- 
quent output  to  the  Nicolet  1180  computer. 

It  would  be  very  interesting  to  excite  TiO  or  ZrO  bands 
with  the  Spectra  Physics  Ar^  laser  and  dye  laser  available 
in  our  laboratory  and  observe  changes  in  the  infrared  spec- 
trum due  to  depopulation  of  the  ground  state.  This  experi- 
ment should  be  particularly  helpful  in  assignment  of  bands 
in  matrices  containing  more  than  one  species.  The  Nicolet 
FT-IR  in  our  laboratory  is  an  excellent  IR  detector  for 
these  experiments,  and  a mobile  matrix  isolation  device  is 
currently  under  construction  for  this  purpose. 

Since  the  lowest  temperature  obtainable  in  this  study 
was  14 °K,  it  was  impossible  to  obtain  Ne  matrices.  The 
construction  of  a liquid  He  cryostat  for  deposition  at  4.2°K 
would  be  an  excellent  addition.  TiO  spectra  in  Ne  matrices 
would  help  resolve  the  problems  of  differentiation  between 
site  effects  and  the  presence  of  other  species  in  the 


matrix . 
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Currently,  an  effort  is  being  made  in  this  laboratory 
to  observe  the  MCD  of  Group  IVb  atoms  in  Ar  matrices.  The 
ZFA  spectra  of  a large  number  of  atomic  species  in  matrices 
have  been  examined  for  a number  of  years,  and  the  matrix 
deposition  conditions  for  many  atomic  species  are  well 
established.  It  would  also  be  interesting  to  look  at  the 
vaporization  products  of  solid  TiO,  ZrO  and  HfO,  which  were 
made  in  this  study  by  vaporization  of  the  corresponding 
dioxide.  The  monoxides  can  be  made  in  bulk  quantities  with 
a purity  exceeding  90%  by  sintering  a mixture  of  the  dioxide 
with  the  element  and  heating  for  several  hours.  The  vapori- 
zation products  of  the  monoxide  obtained  in  this  way  have 
not  been  examined. 


APPENDIX 


The  various  computer  programs  used  in  this  research  are 
presented  below  along  with  brief  comments  on  each  program. 
The  programs  for  the  PET  were  written  in  the  PET  BASIC 
language.  There  were  some  major  drawbacks  to  this  language 
due  to  the  microprocessor  design.  Data  sets  cannot  be 
longer  than  256  points,  and  the  process  of  outputing  data 
to  the  tape  cassette  must  take  into  account  the  operation 
of  the  tape  buffer.  The  programs  for  the  Nicolet  computer 
were  written  in  the  Nicolet  1180  BASIC  language. 

The  first  program,  called  MCDll,  controlled  both 
absorption  and  MCD  experiments,  as  well  as  data  collection 
and  storage.  The  command  RUN30  was  used  to  start  the 
program,  which  is  listed  below. 

MCDll 


10 

J=0 

11 

J=J11 

12 

IF  J<A  GOTO  11 

13 

RETURN 

20 

SAVE" 

MCDll" 

30 

POKE 

4 . 0E4 , 0 

31 

POKE 

3.  3E4,12 

32 

POKE 

4.0E4,4 

33 

POKE 

3. 9E4 , 0 

34 

POKE 

4.9E4,0 

35 

POKE 

4.8E4,255 

36 

POKE 

4.9E4,4 

37 

POKE 

3 . 8E4 , 0 

38 

POKE 

3 . 7E4 , 0 

39 

POKE 

3.8E4,4 
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40  GOTO  350 
50  M=N/4 

60  DIM  A% (250) 

61  DIM  B% (250) 

62  DIM  C% (250) 

63  DIM  D% (250) 

64  DIM  E%  (250) 

65  DIM  F% (250) 

66  DIM  G% (250) 

67  DIM  H%  (250) 

70  FOR  L=2  TO  M 

80  NP=INT(  (L-D/250) 

90  I=L  NP*250 
95  AP=NP+1 

100  ON  AP  GOTO  110,120,130,140,1420144,146,148 
110  A% (I)=PEEK(3. 7E4) 

115  GOTO  150 

120  B% (I)=PEEK(3.7E4) 

125  GOTO  150 

130  C% (I)=PEEK(3.7E4) 

135  GOTO  150 

140  D% (I) =PEEK (3. 7E4) 

141  GOTO  150 

142  E% (I)=PEEK(3.7E4) 

143  GOTO  150 

144  F% (I) =PEEK (3. 7E4) 

145  GOTO  150 

146  G% (I)=PEEK(3.7E4) 

147  GOTO  150 

148  H% (I)=PEEK(3.7E4) 

149  GOTO  150 

150  POKE  3.9E4,8 
160  WL=WL-.25 
170  GOSUB  320 
180  POKE  3.9E4,9 
190  GOSUB  10 

200  FOR  S=1  TO  3 
210  T=0 
220  T=T+1 

230  IF  T<2  GOTO  220 
240  POKE  3.9E4,8 
250  WL=WL-.25 
260  GOSUB  320 
270  POKE  3.9E4,9 
280  GOSUB  10 
290  NEXT  S 
300  NEXT  L 
310  RETURN 

320  V=. 0010497* (WL/2) * (255/7 . 94) +. 5 
325  D==^INT(V) 

330  POKE  4.8E4,D 
335  RETURN 

350  PRINT"DELAY  FACTOR=" ; 

360  INPUT  A 
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370  PRINT" INITIAL  ANGSTROMS="; 

380  INPUT  X 
390  WL=X 
400  GOSUB  320 

410  PRINT"FINAL  ANGSTROMS=" ; 

420  INPUT  Y 
430  N= (X-Y) *4 

440  PRINT" DATA?  YES=1,  NO=0" 

450  INPUT  DC 

460  IF  DC=1  THEN  GOSUB  50 

470  IF  DC=1  GOTO  550 

480  POKE  3.9E4,8 

490  WL=WL-.25 

500  GOSUB  320 

510  POKE  3.9E4,9 

520  N=N-1 

530  GOSUB 

540  IF  N>0  GOTO  480 
550  PRINT"END  OF  SCAN" 

560  PRINT" SAVE  DATA?  YES=1  NO=0" 

600  INPUT  DS 

700  IF  DS=1  THEN  GOSUB  800 
750  END 

800  PRINT" LOAD  DATA  TAPE" 

900  PRINT" INPUT  FILE  NAME" 

950  INPUT  A6$ 

1000  PRINT" INPUT  # DATABLOCKS" 

1100  INPUT  N 

1150  OPEN  1,1,1,A6$ 

1200  FOR  L=1  TO  N 

1300  ON  L GOTO  1400,2000,2700,3400,3910,3920,3930, 
3940 

1400  FOR  I-ITO  250 
1500  X%=A%(I) 

1600  T6=T1 
1700  PRINT#1,X% 

1800  IF  Tl-T6>120  THEN  GOSUB  61900 

1900  NEXT  I 

1950  GOTO  4000 

2000  FOR  J=1  TO  250 

2100  XX=B%(J) 

2200  T6=Tl 
2300  PRINT#1,X% 

2400  IF  Tl-T6>120  THEN  GOSUB  61900 

2500  NEXT  J 

2600  GOTO  4000 

2700  FOR  K=1  TO  250 

2800  X%=C%(K) 

2900  T6=Tl 
3000  PRINT#1,X% 

3100  IF  Tl-T6>120  THEN  GOSUB  61300 

3200  NEXT  K 

3300  GOTO  4000 

3400  FOR  1=1  TO  250 

3500  X%=D% (I) 


165 


3600  T6=T1 
3700  PRINT#1,X% 

3800  IF  Tl-T6>120  THEN  GOSUB  61900 
3900  NEXT  I 
3905  GOTO  4000 

3910  FOR  1=1  TO  250 

3911  X%=E% (I) 

3912  T6=TI 

3913  PRINT#1,X% 

3914  IF  Tl-T6>120  THEN  GOSUB  61900 

3915  NEXT  I 

3916  GOTO  4000 

3920  FOR  1=1  to  250 

3921  X%=F% (I) 

3922  T6=T1 

3923  PRINT#1,X% 

3924  IF  Tl-T6>120  THEN  GOSUB  61900 

3925  NEXT  I 

3926  GOTO  4000 
3930  FOR  1=1  to  250 

3941  X%=G% (I) 

3932  T6=T1 

3933  PRINT#1,X% 

3934  IF  Tl-T6>120  THEN  GOSUB  61900 

3935  NEXT  I 

3936  GO  TO  4000 
3940  FOR  1=1  TO  250 
3940  X%=H%(I) 

3942  T6=T1 

3943  PRINT#1,X% 

3944  IF  Tl-T6>120  THEN  GOSUB  61900 

3945  NEXT  I 
4000  NEXT  L 
4100  CLOSE  1 

4200  PRINT" RIGHT  COMPLETE" 

4300  RETURN 

61900  POKE  59411, 53 :T6=T1 

61910  IF  Tl-T6<6  THEN  61910 

61920  POKE  59411, 61 zRETURN 

62990  PRINT"TYPE  ANY  KEY  WHEN  READY" 

63010  GET  B6$:IF  B6$=""THEN  63010 

63020  RETURN 


The  second  program  was  called  READ6  and  was  used  to 
retrieve  both  absorption  and  MCD  data  from  tape  and  to 
output  this  data  to  the  Nicolet  1180  computer.  READ6  is 
listed  below  and  is  executed  by  the  command  RUN5000. 
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10 

11 

12 

13 

5000 

5001 

5002 

5003 

5004 

5005 

5006 

5007 
5100 
5200 
5300 
5400 
5500 
5600 
5700 
5800 
5900 

6000 

6100 

6200 

6250 

6300 

6400 

6500 

6550 

6600 

6700 

6800 

6900 

7000 

7100 

7200 

7205 

7210 

7215 

7220 

7230 

7240 

7245 

7250 

7255 

7260 

7265 

7270 

7275 

7280 


READ  6 


J=0 

J=J+1  Delay  Subroutine 

IF  J<A  GOTO  11 

RETURN 

DIM  A% (250) 

DIM  B% (250) 

DIM  C% (250) 

DIM  D% (250) 

DIM  E% (250) 

DIM  F% (250) 

DIM  G% (250) 

DIM  H%  (250) 

PRINT" LOAD  DATA  TAPE" 

PRINT" INPUT  FILE  NAME" 

INPUT  A6$ 

PRINT" INPUT  FILE  NAME" 

PRINT" UP  TO  8 BLOCKS  TO  BE  READ" 

INPUT  N 

OPEN  1,1,0,A6$ 

ON  L=1  TO  N 

ON  L GOTO  6000,6300,6600,7000,7210,7240,7260, 
7280 

FOR  1=1  TO  250 
INPUT#1,A% (I) 

NEXT  I 
GOTO  7300 
FOR  J=1  TO  250 
INPUT#1,B% (J) 

NEXT  J 
GOTO  7300 
FOR  K=1  TO  250 
INPUT#1,C% (K) 

NEXT  K 
GOTO  7300 
FOR  1=1  TO  250 
INPUT#1,D% (I) 

NEXT  I 
GOTO  7300 
FOR  1=1  TO  250 
INPUT#1,E% (I) 

NEXT  I 
GOTO  7300 
FOR  1=1  TO  250 
INPUT#1,F% (I) 

NEXT  I 
GO  TO  7300 
FOR  1=1  TO  250 
INPUT#1,G%  (I) 

NEXT  I 

GOTO  7300 

FOR  1=1  TO  250 
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7285  INPUT#1,H% (I) 

7290  NEXT  I 
7300  NEXT  L 

7400  PRINT"END  OF  DATA  READ" 

7500  POKE  4.0E4,0 
7600  POKE  3.9E4,127 
7700  POKE  4.0E4,4 
7800  POKE  3.9E4,0 
7850  POKE  4.9E4,0 
7860  POKE  4.8E4,255 
7870  POKE  4.9E4,4 

8000  PRINT" INPUT  # POINTS  TO  BE  OUTPUT" 
8100  INPUT  NUM 

8110  PRINT" INPUT  DELAY  FACTOR  WHEN  1180" 

8111  PRINT" SAYS  TO  OUPUT  DATA  FROM  PET" 
8120  INPUT  A 

8150  POKE4 . 8E4 , A% (I) 

8160  POKE3.9E4,8 
8200  FOR  L=1  TO  NUM 
8300  NP=INT(  (L-D/250) 

8400  I=L-NP*250 
8800  GO  TO  9500 


8900 

POKE 

4 . 8E4 ,B% 

(I) 

9000 

GOTO 

9500 

9100 

POKE 

4 . 8E4 ,C% 

(I) 

9200 

GOTO 

9500 

9300 

POKE 

4 . 8E4 , D% (I) 

9400 

GOTO 

9500 

9410 

POKE 

4 . 8E4 ,E% 

(I) 

9420 

GOTO 

9500 

9430 

POKE 

4 . 8E4 ,F% 

(I) 

9440 

GOTO 

9500 

9450 

POKE 

4 . 8E4 ,G% 

(I) 

9460 

GOTO 

9500 

9470 

POKE 

4 . 8E4 ,H% (I) 

9480 

GOTO 

9500 

9500 

POKE 

3.9E4,9 

9650 

POKE 

3.9E4,8 

9700 

NEXT 

L 

11700 

END 

The  program,  PEEK,  was  used  to  input  data  from  the  PET 
into  the  Nicolet  1180  using  the  A/D  converter  in  the  Lab  11 
acquisition  mode.  The  program  stored  the  data  on  the 
Diablo  disk  system  using  an  appropriate  disk  directory 
filename.  The  incoming  data  was  displayed  on  the  digital 
oscilloscope,  and  the  data  could  also  be  plotted.  The  data 
input  program,  PEEK,  is  listed  below. 
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PEEK 


100  DIM  X(2001) ,Y(2001) 

110  PRINT"WANT  TO  INPIT  FROM  PET  NO-0" 

111  INPUTQO 

112  IFQ0=0 

113  CALL  FCLEAR(X, 2001) 

114  CALL  FCLEAR(Y, 2001) 

150  PRINT" INPUT  #POINTS  TO  ACQUIRE" 

160  INPUT  N5 
162  LET  N4=N5*2+2 

165  PRINT  "OUTPUT,  DATA  FROM  PET" 

170  CALL  FAQT(Y,N5, 200, 200, 1036, 0,1) 

172  PRINT"ACQUIRE  DONE" 

175  PRINT  "INPUT  DATA  FILENAME" 

176  CALL  INFILE 

177  CALL  FSTORE(Y,N4) 

180  PRINT  " WANT  TO  PLOT?  NO=0" 

181  INPUT  Q1 

182  IFQ1=0  THEN  110 

200  PRINT"INPUT  NO  POINTS,  X AXIS  LEN,Y  AXIS  LEN" 

300  INPUT  N,L0,L1 

700  CALL  DISPOF 

800  CALL  FDISPI (Y, 2000, 200) 

1500  CALL  PLOTS 

1501  PRINT" INPUT  PLOT  ORIGIN :X,Y" 

1502  INPUT  X9,Y9 

1503  LET  Y9=Y9+0.5 

1510  CALL  FDPLOT  (X9,Y9,-3) 

1520  PRINT" INPUT  STARTING  WVELENGTH" 

1530  INPUT  W1 
1535  LET  Nl=-1 
1540  FOR  I=0TO  Nl 
1550  LET  X(I)=W1-I 
1561  NEXT  I 
1570  GOSUB  15000 
1600  LET  X1$="ANG" 

1610  LET  Y1$="ABS" 

1614  PRINT  "7VFTER  FSCALE" 

1615  PRINT  "X(N) ,X(N+1)=";X(N) ,X(N+1) 

1620  CALL  FAXIS(0,0,X1$,3,L0,0,X(N) ,X(N+1) ,3,3) 

1630  CALL  FAXIS (0, 0, Yl$ , 3,L1, 90, Y (N) , Y (N+1) , 3, 3) 

1640  CALL  FLINE(X,Y,N, 1,0,0) 

1650  PRINT" PLOT  DONE" 

1655  PRINT"WOULD  YOU  LIKE  TO  PLOT  IT  DIFFERENTLY  NO=0)" 

1656  INPUT  Q2 

1657  IF  Q2=0  THEN  110 

1658  GO  TO  200 
1700  GO  TO  110 

15000  LETX(N)=1E30 

15001  LET  Y(N)=1E30 
15010  LET  M1=-1E30 
15020  LET  M2=-1E30 
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15030  FOR  I=OTO  N-1 

15040  IF  X(I)>X(N)  THEN  15060 

15050  LET  X(N)==X(I) 

15060  IFX(I)<M1  THEN  15080 
15070  LET  M1=X(I) 

15080  IFY(I)>Y(N)  THEN  15100 
15090  LETY(N)=Y(I) 

15100  IF  Y(I)<M2  THEN  15120 
15110  LET  M2=Y(I) 

15120  NEXT  I 

15130  LET  X(N+1)=(M1-X(N) )/L0 
15140  LET  Y(N+1)=(M2-Y(N) )/Ll 
15150  RETURN 
20000  END 


The  program,  SPECAS,  was  used  to  assemble  data  sets 
stored  on  the  disk  into  full  spectra.  Calibration  and 
wavelength  overlap  of  data  sets  was  performed  as  well  as 
storage  and  plotting  of  the  final  spectrum.  This  spectrum 
assembler  program,  SPECAS  is  listed  below. 


SPECAS 


100  DIM  Y1 (2002)  ,Y2  (1000)  ,Y(6002 

150  CALL  FCLEAR(Y, 6002) 

151  CALL  FCLEAR(Y1,2002) 

152  CALL  FCLEAR(Y2,1000) 

200  PRINT" INPUT  NUMBER  OF  POINTS  IN  DATA  BLOCK" 

300  INPUT  N 

301  LET  Nl=(N*2)+2 

310  PRINT" INPUT  DATA  FILE  NAME  OF  INCOMING  DATA" 

320  CALL  INFILE 

330  CALL  FL0AD(Y1,N1) 

331  PRINT  "IS  THIS  AN  ASSEMBLY  OF  ABSORPTION  OR  MCD 
RUNS?  ABS=A  MCD=M" 

332  INPUT  BIS 

333  IF  B1$="A"  THEN  2050 

335  FOR  1=0  TO  N-1 

336  LET  Y1 (I)=Y1(I)-1535 

337  NEXT  I 

340  CALL  DISPOF 

350  CALL  FDISPI (Y1,N, 200) 

400  PRINT" INPUT  ZERO  FIELD" 

450  INPUTZ2 

455  LET  Z2=Z2-1535 

460  PRINT" INPUT  CALIBRATION  FACTOR" 

470  INPUT  C2 

475  LET  C2=C2-1535 
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1200  PRINT"INPUT  SENSITIVITY  FOR  STANDARD" 

1300  INPUT  SI 

1400  PRINT"INPUT  SENSITIVITY  FOR  DATA" 

1500  INPUT  D1 

1600  PRINT" INPUT  DELTA  ABSORBANCE  FOR  STANDARD" 

1700  INPUT  S2 

1800  LET  Fl= (Dl*S2)/ (S1*C2-Z2) 

1900  FOR  1=0  TO  N-1 

1950  LET  Y1(I)=(Y1(I)-Z2) *F1 

2000  NEXT  I 

2050  CALL  DISPOF 

2060  CALL  FDISPI (Y1,N,200) 

2100  PRINT" INPUT  STARTING  WAVELENGTH" 

2300  INPUT  W1 

2310  PRINT"WANT  OVERLAP  ADJUSTMENT  NO=0 

2311  INPUT  Q1 

2312  IF  Q1=0  THEN  2395 

2313  FOR  1=0  TO  5999 

2314  LET  Y(I)=Y(I)/1E5 

2315  NEXT  I 

2320  PRINT" INPUT  SHORTEST  WAVELENGTH  IN  OVERLAPPING 
REGION" 

2330  INPUT  W2 
2335  LET  S1=0 

2339  LET  W3=Wl-W2 

2340  FOR  1=0  TO  W3 

2350  LET  S1=S1+Y (8500-Wl+I) -Y1 (I) 

2360  NEXT  I 

2370  LET  A1=S1/(W3+1) 

2371  PRINT"AVG  DIF.  IN  OVERLAP=";Al 
2380  FOR  1=0  TO  W3 

2384  LET  Y2 (I) =Yl (I) 

2385  LET  Y1 { I ) =Yl ( I ) +Al 
2390  NEXT  I 

2395  LET  N1=N-1 
2400  FOR  1=0  TO  Nl 
2500  LET  J=8500-W1+I 
2600  LET  Y (J)  =Y1 (I) 

2700  NEXT  I 

2701  FOR  I-O  to  5999 

2702  LET  Y(I)  =Y(I) *1E5 
2303  NEXT  I 

2710  CALL  DISPOF 

2720  CALL  FDISPI (Y , 6 002 , 200 ) 

2725  PRINT"  DOES  OVERLAP  LOOK  OKAY?  YES=Y" 

2730  INPUT  Q$ 

2735  IFQ$="Y"  THEN  2800 
2740  FOR  1=0  TO  W3 
2750  LET  Y1 (I) =Y2 (I) 

2760  NEXT  I 

2770  PRINT"INPUT  A DIFFERENT  AVG  OVERLAP  DIFFERENCE" 

2780  INPUT  A1 

2790  GO  TO  2380 

2800  PRINT"MORE  DATA  NO=N" 
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700  PRINT" INPUT  STARTING 
800  INPUT  W1 

805  PRINT"WOULD  YOU  LIKE 

806  INPUT  Ql$ 

807  IF  Q1$="Y"  THEN  810 

808  GO  TO  900 

810  PRINT  "INPUT  LONGEST 
820  INPUT  W3,W2 
825  FOR  1=0  TO  W3-W2 
830  LET  Y(I)=Y(W1-W3+I) 
840  NEXT  I 
850  LET  N=W3-W2+1 
860  LET  W1=W3 
900  FOR  1=0  TO  N-1 
1000  LET  X(I)=W1-I 
1100  NEXT  I 
1200  PRINT"LIKE  SAME 
1210  INPUT  Q2$ 

1220  IF  Q2$="N"  THEN 
1230  LET  X(N)=X3 
1235  LET  X(N+1)=X4 
1240 
1245 
1250 
1300 
1400 
1450 
1500 
1600 
1601 
1602 
1603 
1610 
1615 
1620 
1630 

1700 

1701 

1702 

1703 

1705 

1706 

1707 
1710 
1720 
1730 
1746 
1750 
1760 
1770 
1780 
1800 

15000 

15001 


WAVELENGTH" 

TO  PLOT  A SMALL  PART?  YES=Y 
AND  SHORTEST  WAVELENGTHS" 


SCALE  FACTOR  AS  LAST  PLOT?  YES=Y 
1300 


LET  Y(N)=Y3 
LET  Y(N+1)=Y4 
GO  TO  1700 
LET  X1$="CM-1" 

LET  Y1$="ABS" 

LET  X2$="ANG" 

PRINT" INPUT  X,Y  AXIS  LENGTHS" 

INPUT  L0,L1 

PRINT"INPUT  ORIGN  X,Y" 

INPUT  01,02 

CALL  FDPLOT(01,02,-3) 

GOSUB  15000 
LET  L3=0 

CALL  FAXIS(0,L3,X2$,3,L0,0,X(N) ,X(N+1) ,4,4) 
CALL  FAXIS(0,0,Y1$,3,L1,90,Y(N) ,Y(N+1) ,4,2) 
CALL  FLINE(X,Y,N, 1,0,0) 

PRINT"WANT  WAVENUMBER  SCALE?  NO=N" 

INPUT  T2$ 

IF  T2$="N"  THEN  20000 

FOR  1=0  TO  N-1 

LET  X(I)=l/( (Wl-I) *lE-8) 

NEXT  I 

LET  R=0 

GOSUB  15000 

LET  Wl=l/ (Ml*lE-8) 

LET  W=-X(N+1) 

CALL  FAXIS ( 0 , 0 , XI $ , 4 , LO , 0 , W1 , W , 4 , 3 ) 

CALL  FDPLOT(0,L3, 3) 

CALL  FDPLOT(0,0,2) 

CALL  FDPLOT (LO, 0, 2) 

GO  TO  20000 
LET  X(N)=1E30 
LET  Y(N)=1E30 
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15010 

15020 

15030 

15040 

15050 

15060 

15080 

15090 

15100 

15110 

15120 

15121 

15122 

15123 

15125 

15126 

15127 
15130 
15140 
15150 
15160 
15170 
15180 
15190 
15200 
20000 
20001 
20002 
20005 


LET  Ml=-1E30 
LET  M2=:-1E30 
FOR  1=0  TO  N-1 
IF  X(I)>X(N)  THEN  15060 
LET  X(N)=X(I) 

IF  X(I)<M1  THEN  15080 
IF  Y (I) >Y (N)  THEN  15100 
LET  Y(N)=Y(I) 

IF  Y(I)<M2  THEN  15120 
LET  M2=Y(I) 

NEXT  I 

IF  R=0  THEN  15130 

PRINT" INPUT  TRUNCATION  FACTOR" 

INPUT  Tl 

LET  X9=X(N)/T1 

LET  X(N)=T1*INT(X9) 

LET  M1=INT(M1/T1) *T1+T1 
LET  X(N+1)  = (M1-X(N)  )/L0 
LET  Y(N+1)  = (M2-Y(N)  )/Ll 
LET  X3=X(N) 

LET  X4=X(N+1) 

LET  Y3=Y(N) 

LET  Y4=Y(N+1) 

LET  R=1 
RETURN 

PRINT"W7^T  TO  PLOT  MORE  DATA?  YES=Y" 
INPUT  Sl$ 

IF  S1$="Y"  THEN  60 
END 
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